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I.  INTRODUCTION 

This  report  describes  the  development  of  a  model  of  soil  breakdown  based 
on  a  large  number  of  experiments  which  were  performed  under  the  auspices  of 
the  Air  Force  Weapons  Laboratory  (AFWL). 

The  purpose  of  this  report  is  to  make  use  of  this  experimental  data  to 
develop  a  phenomenological  model  of  soil  breakdown  suitable  for  use  in  a  com¬ 
puter  code.  The  experiments  performed  by  Mission  Research  Corporation  (MRC) 
have  been  used  because  of  the  detail  of  spatial  voltage  variation  in  the  pres¬ 
ence  of  soil  breakdown  that  was  obtained.  Nonpurturbing  voltage  probes  were 
inserted  into  the  soil  samples  to  measure  voltage  versus  time  at  a  large 
number  of  discrete  points. 

Soil  breakdown  characterized  by  corona  formation  and  streamer  propagation 
has  been  recognized  for  several  decades  (Ref.  1).  It  was  assumed  that  a  corona 
developing  about  a  buried  electrode  (subjected  to  a  large  pulse)  would  have 
streamers  forming  within  it  (Ref.  2).  Corona  is  discussed  in  Reference  3  and 
4.  An  excellent  reference  on  both  linear  and  nonlinear  soil  properties  from 
this  era  is  found  in  Reference  5.  In  Reference  6,  B.  E.  Lewis  has  tried  to 
model  the  decay  of  the  corona  as  the  pulser  voltage  relaxes.  If  the  models 
of  this  early  work  are  somewhat  modified,  they  will  fit  the  experimental  data 
now  in  use.  These  models  have  been  quantified  in  such  a  way  that  they  can  be 
uded  in  computer  codes. 

It  may  be  useful  to  read  Appendix  A  at  this  point  before  tackling  the 
material  of  this  report.  This  appendix  is  a  description  of  the  nature  of 
soil  from  a  geologist's  point  of  view.  This  description  defines  the  com¬ 
ponents  of  soil  and  then  decribes  through  Archie's  law  how  each  component 
contributes  to  the  linear  ohmic  conductivity  of  soil.  It  also  provides  a  good 
basis  on  which  to  build  an  understanding  of  the  nonlinear  behavior  of  soil. 

In  Section  II,  the  relevant  feature  of  the  experimental  configurations 
which  was  used  In  the  experiments  are  discussed.  Section  III  is  an  overview 
of  the  model;  Section  IV  is  a  development  of  the  electrostatic  model  while 
Section  V  develops  the  dynamic  model  of  relaxation.  Section  VI  discusses 
streamers  and  Section  VII  is  a  summary  of  the  important  points  made  in  this 
report. 
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II.  EXPERIMENTAL  SYSTEMS 


Nearly  all  this  report  is  based  upon  experiments  which  were  performed  by 
MRC.  This  work  was  performed  in  two  experimental  configurations  both  of  which 
contained  cylindrical  symmetry: 

a.  A  cylindrical  wedge  of  angle  it/2.  The  thickness  of  the  wedge 
was  either  0.1  m  or  0.3  m  and  the  radius  was  0.5  or  1.0  m.  When  an  experi¬ 
ment  is  discussed,  the  wedge  thickness  and  radius  will  be  stated.  Experiments 
will  henceforth  be  called  Experiments  A  or  Exp.  A. 

b.  A  full  2  ir  cylinder.  The  length  was  1  m  and  the  radius  0.5  m 
was  1  m.  Experiments  performed  in  this  configuration  will  be  called  Exp.  B. 

The  MRC  performed  the  Exp.  A  at  Physics  International,  San  Leandro, 
California.  The  basic  configuration  is  shown  in  Figure  1.  The  soil  sample  is 
in  the  shape  of  a  ir/2  wedge  with  a  center  conductor  of  radius  a,  an  outer  con¬ 
ductor  (a  copper  sheet)  of  radius  b,  and  a  thickness  L.  In  some  of  the 
experim^hts,  the  outer  conductor  was  connected  to  pulser  ground  through  a 
resistor  of  200  a.  The  top  and  bottom  of  the  soil  sample  as  well  as  the  sides 
at  0  and  w/2  rad  were  of  an  insulating  material.  The  intervening  space  was 
filled  with  soil.  For  the  experiments  of  interest,  a  =  0.01  m,  b  =  0.5,  1  m, 
and  L  =  0.1,  0.3  m. 

The  center  conductor  was  driven  by  a  30  to  100  kV,  0.9  yF  pulser.  The 
interesting  feature  about  these  experiments  was  that  the  voltage  in  the  soil 
sample  was  measured  at  discrete  points  as  a  function  of  time  by  an  array  of 
voltage  probes.  A  probe  had  a  high  impedance  between  its  small  spherical  tip 
(in  contact  with  soil)  and  ground.  Hence,  it  drew  little  current.  The  probes 
were  designed  so  that  they  would  not  purturb  significantly  the  distribution  of 
voltage  in  the  soil  sample. 

The  array  of  probe  positions  which  were  available  are  shown  in  Figure  1. 
There  were  10  data  channels  available;  nine  of  these  were  used  for  monitoring 
probe  voltages  while  the  tenth  was  used  to  monitor  total  current  I  through  the 
sample.  Thus,  for  a  given  discharge  of  the  pulser  (shot),  nine  probes  could 
be  monitored.  The  probes  were  either  located  in  the  nine  radial  positions  in 
the  figure  or  in  the  nine  azimuthal  positions. 
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The  soil  types  used  in  Exp.  A  which  are  relevant  are: 

a.  Sand  with  an  ohmic  resistivity  pQ  of  about  150  to  2000  n-m; 

b.  a  sand  and  gravel  composite  labeled  SOIL-B;  a  typical  pQ  is  130  fi-m. 

The  ffllC  performed  Exp.  B  at  the  McDonnel 1 -Douglas  Lightning  Facility 
(MCAIR),  St.  Louis,  Missouri.  The  sample  was  a  full  2  it  cylinder  with  a  central 
conducting  rod  of  radius  a.  The  outer  conductor  was  a  standard  steel  culvert 
of  radius  b.  The  length  of  the  system  was  L.  The  system  had  a  top  and  bottom 
plate  which  were  of  an  insulating  material  and  the  intervening  space  was  filled 
with  soil.  For  the  experiments  in  this  report;  a  =  0.01  m,  b  =  0.5  m  and 
L  =  1  m. 


The  center  conductor  was  driven  by  one  or  two  five-stage  pulsers  at  the 
MCAIR  facility.  Each  Marx  pulser  could  be  charged  to  15  to  24  kV  per  stage 
with  a  total  capacitance  of  230  uF  divided  between  the  number  of  stages  used. 
The  same  high  impedance  voltage  probes  were  used  in  these  experiments  that 
were  used  in  Exp.  A.  The  probes  were  inserted  into  the  soil  through  the  top 
and  bottom  insulating  plates  in  either  a  radial  or  azimuthal  array.  Of 
interest  here  is  an  azimuthal  array  described  in  more  detail  in  Section  VI. 
Twenty-six  data  channels  were  available.  The  soil  of  interest  in  this  report 
is  300  fl-m  sand. 

Finally,  there  is  a  set  of  experiments  called  Exp.  C.  Which  were  per¬ 
formed  by  JAYC0R.  These  were  point-plane  experiments  where  the  soil  resis¬ 
tivity  was  50  fi-m. 
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III.  THE  MODEL,  AN  OVERVIEW 


The  model  that  is  adapted  in  the  report  consists  of  two  components:  the 
corona  and  streamers.  Although  these  components  are  inseparable,  in  pre¬ 
senting  an  overview  of  the  model,  they  will  be  discussed  individually. 

1 .  CORONA 

In  discussing  the  corona,  Exp.  A  will  be  used.  In  Exp.  A,  the  voltage 
V(r)  with  respect  to  ground  was  measured  as  a  function  of  time  by  high  imped¬ 
ance  probes  at  the  radial  locations  shown  in  the  cylindrical  wedge  geometry 
illustrated  in  Figure  1.  Also  measured  was  the  total  current  I  flowing 
through  the  wedge  as  a  function  of  time.  Thus,  it  is  possible  to  plot  the 
ratio  V/I  (the  effective  resistance  between  the  location  of  the  probe  and 
ground)  as  a  function  of  radial  position  r,  or  more  conveniently,  log  r.  For 
a  typical  shot,  this  plot  is  shown  in  Figure  2. 

If  the  cylindrical  symmetry  is  assumed,  the  cylindrical  wedge  is  electri¬ 
cally  divided  into  two  regions  (Fig.  2).  The  division  between  the  two  regions 
occurs  at  a  radius  rQ,  which  for  the  parameters  of  Figure  2  would  be  about 
0.16  m.  In  each  region,  V/I  and,  hence,  V  is  nearly  proportional  to  log  r; 
since  in  this  symmetrical  case,  E  =  r  Er,  we  are  then  led  to  Er  being  propor¬ 
tional  to  r  1  in  each  region,  although  the  proportionality  constants  are  dif¬ 
ferent.  If  L  is  the  thickness  of  the  wedge,  then  the  current  density  is 


where  p  is  the  soil  resistivity  and  A  is  the  proportionality  constant  for  a 
region.  Equation  1  states  that  p  is  constant  in  each  region  and  the  dif¬ 
ference  in  proportionality  constants  can  be  accounted  for  by  the  resistivity 
being  different  in  each  region. 

If,  as  in  Figure  1,  a  is  the  radius  of  the  central  high  voltage  conductor 
and  b  is  the  radius  of  the  outer  grounded  cylinder,  then  the  following  is 
assumed : 


SHOT0399 


Figure  2.  Resistance  (kn)  for  shot  399  at  various  times.  Curve  A,  second  from  top;  B,  third;  C,  fifth 
D,  sixth;  E,  fourth;  F,  top  curve.  Total  current  through  sample  shown  at  bottom. 
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a.  For  rQ  £  r  _<  b,  the  soil  is  not  broken  and  is  an  ohmic  state 
such  that  p  =  p0,  the  ohmic  value;  furthermore,  at  r  =  rQ,  Ep  =  Eg,  a  critical 
electric  field  above  which  the  soil  will  break  down  but  below  which  it  will 
not;  this  is  called  the  ohmic  region. 

b.  For  a  £  r  £  rQ,  Er  >  Eg  and  the  soil  is  broken  down  in  some 
fashion;  here,  p  =  px  <  pQ.  This  is  called  the  corona  region. 

This  two  region  description  of  soil  breakdown  is  referred  to  as  the 
electrostatic  model.  A  soil  sample  which  is  electrically  stressed  such  that 
the  electrostatic  model  applies  is  said  to  be  in  the  electrostatic  phase. 

Although  it  has  been  assumed  that  the  corona  will  adapt  the  symmetry  of 
the  electrode  system  e.g.,  cylind  „al  or  spherical,  assymetries  in  the  soil 
or  other  perturbations,  e.g.,  st  mers,  may  cause  the  corona  to  depart  from 
the  symmetry  of  the  electrode  sy  m.  The  electrostatic  model  may  be  general¬ 
ized  to  be  that  in  which  the  ele  .<*  illy  stressed  soil  is  still  divided  into 
two  regions  as  illustrated  in  Figure  3.  In  the  region  containing  the  high 
voltage  electrode  or  the  corona  region,  the  resistivity  is  px  (r)  <  pQ.  In 
the  region  containing  the  grounded  electrode  or  the  ohmic  region,  the  soil  is 
not  broken  down  and  the  resistivity  is  just  pQ.  Furthermore,  on  the  boundary 
between  the  two  regions,  in  the  ohmic  region,  the  electric  field  is  not  neces¬ 
sarily  normal  to  the  boundary,  but  its  magnitude  is  constant  everywhere  on 
the  boundary  and  is  Eg,  the  critical  electric  field.  Finally,  a  critical 
current  density  can  be  defined  at  the  boundary  in  the  ohmic  region, 

Oc  -  Vpo  (2) 

The  magnitude  of  this  current  density  is  constant  everywhere  on  the  boundary. 

It  is  of  interest  to  see  what  precedes  the  electrostatic  phase  and  how 
this  phase  decays  as  the  voltage  on  the  high  voltage  electrode  decays.  Call 
V0  the  voltage  on  the  central  high  voltage  electrode.  As  V0  increases,  the 
soil  will  remain  ohmic  until  the  electric  field  at  the  electrode  Ep(a)  =  E0 
approaches  about  1  MV/m.  Assuming  a  rise  time  of  V0  of  the  order  of  a  micro¬ 
second,  there  will  be  a  delay  time  t^  until  the  soil  starts  to  break  down  and 
the  soil  remains  ohmic.  During  this  phase,  it  has  been  observed  that  the 
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CO 


Boundary  between  an  corona  region  and  ohmic  soil  region; 
and  J  =  E„/p  . 


resistivity  of  the  soil  actually  is  larger  than  pQ.  For  example,  in  shot  012, 
while  the  ohmic  resistance  of  the  soil  is  14.7  left  determined  by  low  voltage 
measurements,  during  the  prebreakdown  phase,  the  resistance  was  measured  to  be 
16.1  kQ  (Exp.  A). 

Once  breakdown  begins,  the  resistance  of  the  soil  drops  rapidly  in  time 
as  the  total  current  through  the  sample  I  rises  even  though  the  V0  is  dropping. 
After  a  few  tens  of  ys,  rQ  begins  to  appear  as  the  sample  begins  to  enter  the 
electrostatic  phase.  Two  cases  must  be  distinguished  at  this  point.  The 
first  is  that  VQ  can  be  sufficiently  large  to  cause  the  corona  to  extend  all 
the  way  to  the  grounded  conductor  at  r  =  b.  We  call  this  corona  attachment. 
Then  I  will  increase  to  a  large  value,  and  will  only  be  limited  by  the  low 
resistance  of  the  soil.  Corona  attachment  can  occur  if  VQ  is  so  large  that 
for  a  symmetric  corona,  rQ  =  b,  or  rQ  exceeds  a  critical  radius.  Breakdown  to 
the  outer  conductor  can  also  occur  by  streamer  attachment  discussed  in 
paragraph  2.  The  second  case  is  that  if  VQ  is  not  large  enough  for  any  of 
these  events  to  occur,  the  corona  remains  stable,  I  will  increase  to  a  maximum 
value,  and  then  decline  to  zero  as  the  voltage  VQ  falls  to  zero.  For  a  sym¬ 
metric,  stable  corona,  see  in  Figures  B.21  and  B.22  (Appendix  B)  how  the 
current  I  rises  as  the  voltage  V0  falls  after  the  onset  of  breakdown  for  shot 
399. 

For  shot  399  at  30  ys,  the  ohmic  region  could  be  distinguished  as  can  be 
seen  in  Section  IV  and  there  is  indeed  a  region  of  corona  where  px  <  pQ.  The 
current  reaches  its  maximum  about  200  ys.  From  30  to  200  ys,  the  electro¬ 
static  phase  is  clearly  defined.  Thus,  it  takes  the  order  of  30  ys  for  the 
electrostatic  phase  to  be  reached,  which  is  typical  for  shots  of  Exp.  A. 

If  Jc  is  nearly  normal  to  the  boundary  separating  the  corona  from  the 
ohmic  region,  and  since  its  magnitude  is  constant  over  this  surface,  then 


so  that  for  Eg  being  regarded  as  a  constant,  the  surface  of  the  boundary  is  a 
maximum  when  I  is  a  maximum,  for  the  symmetric  corona,  rQ  is  a  maximum  when  I 
is  a  maximum.  Maximum  I  and  rQ  should  depend  on  maximum  V0,  and  as  maximum  or 
peak  VQ  increases,  so  does  maximum  rQ. 
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As  the  current  decays  after  reaching  its  maximum,  then  the  symmetric 
cr-ona  enters  what  is  called  the  dynamic  phase.  Modeling  the  decay  of  the 
corona  after  Liew,  et  al .  (Ref.  6)  a  third  region  is  postulated  to  exist, 
bounded  by  an  inner  radius  rL ,  where  the  current  density  if  Jc,  and  an  outer 
radius  rQ,  corresponding  to  rQ  at  maximum  I. 

The  three  regions  of  the  dynamic  mode  are  defined  as  follows: 

Region  1:  a  <:  r  £  rx ,  J  >_  Jc  with  J  =  Jc  at  r  =  rx ;  p!  =  const. 

Region  2:  ri  <.  r  <.  r0»  J  <.  Jc  with  P2  dependent  on  r. 

Region  3:  rQ  £  r  £b,  the  ohmic  region  with  p  =  pQ. 

The  dynamic  model  is  discussed  in  Section  V. 

2.  STREAMERS 

It  is  well  established  that  a  second  method  of  soil  breakdown  occurs  in 
which  narrow  channels  of  high  conductivity  form  at  the  high  voltage  electrode 
and  proceed  to  grow  outward  at  some  velocity.  These  channels  have  been  called 
streamers.  A  model  proposed  in  Reference  2  has  been  adapted  in  which  a 
streamer  is  viewed  as  contained  inside  the  corona  as  shown  in  Figure  4.  In 
this  model,  the  corona  forms  first,  the  streamer  may  initiate  at  some  point  on 
the  high  voltage  electrode  for  reasons  that  are  not  understood.  The  streamer 
grows  outward  distorting  the  symmetry  of  the  corona.  In  this  figure  a  sym¬ 
metric  corona  (dashed  line)  containing  no  streamers  is  contrasted  with  one 
containing  one  or  more  streamers.  The  streamers  always  stay  inside  the 
corona's  boundary  and  the  electric  field  remains  at  Eg  at  the  boundary  of  the 
corona. 

This  model  is  reasonable  if  one  accepts  that  the  resistivity  of  the 
streamer  is  very  low.  Then  the  voltage  drop  along  the  streamer  will  be  small, 
and  its  potential  will  approximate  that  of  tne  high  voltage  electrode.  Just 
as  this  electrode  will  cause  soil  breakdown  some  distance  out  from  it,  so  the 
streamer,  which  must  have  a  large  electric  field  at  its  surface,  will  cause 
breakdown  in  the  nearby  surrounding  soil.  This  model  of  a  corona  surrounding 
a  streamer  has  been  adopted  by  Longmire  et  al . 

It  is  also  predicted  by  the  University  of  Arizona  that  once  a  streamer  is 
formed,  the  thermodynamics  and  electrical  properties  of  the  streamer  tip  will 
cause  it  to  propagate  outward  with  a  predictable  velocity.  Eventually,  the 
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Figure  4.  A  buried  spherical  electrode  with  streamers  propagating  outward  from  points  A  and 
B.  Streamers  are  embedded  in  corona  (outer  medium  dark  line).  Dashed  semicircle 
is  corona  boundary  if  streamers  are  absent.  Current  into  electrode  increases  from 
(a)  to  (d)  as  does  corona  boundary  surface.  (Fig.  14  of  Ref*  2). 
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streamer  will  push  the  corona  surrounding  its  tip  out  to  the  grounded  electrode 
at  r  =  b,  followed  by  the  streamer  itself  attaching  to  this  electrode.  This 
form  of  breakdown  to  the  grounded  conductor  is  called  streamer  attachment. 

In  Section  VI,  we  are  able  to  make  an  estimate  of  the  radius  of  a 
streamer  that  occurred  in  a  nearly  symmetric  corona.  In  that  case,  it  was 
found  that  the  radius  was  about  1  cm. 
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IV.  THE  ELECTROSTATIC  MODEL  OF  THE  CORONA 

As  seen  in  the  preceding  section,  the  electrostatic  model  is  the  cor¬ 
nerstone  of  the  corona  model.  In  this  section,  some  other  possible  corona 
models  will  be  discussed  and  the  use  of  the  electrostatic  model  justified. 

The  details  of  the  theory  of  the  electrostatic  model  will  be  presented  first 
for  the  symmetric  cylindrical  case  and  then  for  the  symmetrical  spherical 
case.  Next,  the  electrostatic  phase  through  experimental  data  will  be  charac¬ 
terized.  Finally,  it  will  be  shown  that  the  peak  current  can  be  accurately 
predicted  by  the  electrostatic  model. 

1.  OTHER  MODELS 

In  referring  to  Figure  2,  since  V(r)  appears  to  vary  in  a  manner  propor¬ 
tional  to  log  r  in  both  the  corona  and  ohmic  regions,  it  is  concluded  that  in 
each  region  Ep  is  proportional  to  r  Thus,  it  is  assumed  that 

V(r)  -  V0  -  Vx  log  (r/a)  (4) 

On  the  other  hand,  Longmire  el  al .  claim  that  the  electric  field  in  the 
corona  is  clamped  to  that  required  for  electron  avalanching  to  be  just 
balanced  by  electron  diffusion  and  attachment  (Ref.  7).  This  electric  field 
is  estimated  to  be  about  1  MV/m.  In  any  case,  Er  is  assumed  to  be  independent 
of  r,  or  Ep  =  const. 

This  would  lead  to  a  linear  dependence  of  V(r)  on  r, 

V(r)  •  V0  ♦  H  (1  -  L)  (5) 

An  attempt  was  made  to  distinguish  between  these  types  of  voltage 
variation  in  the  corona  by  considering  in  detail  the  data  from  shot  399  of 
Exp.  A.  The  raw  data  or  data  base  from  this  shot  are  shown  in  Figures  B.l 
through  B.10  (Appendix  B).  Figures  5  through  8  show  the  plotted  data  from  the 
radial  probes  at  the  locations  shown  in  Figure  1  versus  log  r  at  two  early 
times,  30  and  60  us  and  at  two  times  when  the  total  current  I  was  nearly  its 
maximum,  i.e.,  100  and  300  us. 
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solid  line)  and  linear  r  (dashed  line)  fit  to  radial  voltage  probe  data 
it  399;  t  =  60  ys. 


Figure  7.  Log  r  (solid  line)  and  linear  r  (dashed  line)  fit  to  radial  voltage  probe  data 
for  shot  399;  t  =  100  ys. 


Consider  Figure  7  for  100  us.  The  probes  in  Figure  1  are  numbered  1 
through  9  and  the  voltage  at  each  probe  is  shown  by  a  heavy  dot  in  Figure  7 
with  the  number  of  the  probe  from  which  it  was  obtained  shown  above  the  point. 
The  experimenters  claim  that  the  random  error  in  a  probe  measurement  had  a 
standard  deviation  of  4  percent  which  is  interpreted  to  mean  a  random  error  of 
±2%.  This  error  is  shown  by  the  error  bars  on  the  data  points. 

A  least  squares  logrithmic  regression  fit  has  been  done  of  Figure  5  on 
points  1,  2,  and  3  and  the  resulting  curve  is  shown  at  the  solid  curve  in 

Figure  7.  These  three  points  were  used  since  these  were  the  only  probes  in 

the  corona  region. 

A  least  squares  linear  regression  fit  of  Figure  6  has  been  done  on  points 
1,  2,  3,  4,  and  5  and  the  resulting  curve  is  shown  as  dashed  in  Figure  7. 
Although  the  error  bars  tend  to  support  the  log  r  dependence,  this  conclusion 
is  somewhat  tentative.  There  were  only  three  probes  in  the  corona  region;  had 
there  been  more,  a  more  certain  conclusion  might  have  been  reached.  Notice 
that  to  distinguish  between  the  two  types  of  voltage  dependence  on  r,  very 
precise  measurements  of  V  are  required.  A  ±2%  error  is  the  outside  limit  of 
what  can  be  tolerated. 

The  same  kind  of  least  squares  fit  of  Figure  5  and  6  have  been  made  on 

the  same  groupings  of  points  in  Figures  5,  6,  and  8  for  30,  60,  and  300  us, 

respectively.  At  30  us,  neither  dependence  seems  correct  based  on  the  error 
bars.  At  60  us,  the  log  r  dependence  gives  a  better  fit  than  the  linear  fit 
as  is  the  case  at  100  us.  At  300  us,  the  log  r  fit  is  nearly  exact.  Thus,  as 
time  proceeds  toward  I  reaching  its  maximum,  the  log  r  fit  gets  progressively 
better. 

In  what  follows,  a  log  r  dependence  will  be  assumed.  From  the  point  of 
view  of  trying  to  develop  a  phenomenological  model  of  the  corona  which  is 
suitable  for  computer  codes,  whether  on  selects  Figure  5  or  6,  certain  impor¬ 
tant  parameters  would  be  predicted  in  either  case.  Voltage  v(r)  would  be 
satisfactorily  predicted  since  the  difference  between  the  solid  and  dashed 
curves  In  the  figures  Is  not  large.  Also,  the  total  current  I  and  current 
density  as  well  as  rQ,  the  boundary  between  the  two  regions,  would  be  pre¬ 
dicted.  What  would  be  lost  would  be  Er  and  the  dependence  of  px  on  r. 
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In  Figures  5  through  8  a  least  squares  logrithmic  regression  fit  on  points 
5  through  9  in  the  ohmic  region  has  been  performed.  The  fit  is  excellent. 

The  slope  of  this  curve  determines  the  ohmic  resistivity.  It  can  be  seen 
that  there  is  a  significant  apparent  variation  of  p0  among  the  four  times. 
Assume  that  this  variation  is  not  real,  but  reflects  errors  in  measuring  I, 
which  was  used  in  the  computation  of  pQ. 

The  intersection  of  the  log  r  curve  for  points  in  the  ohmic  region  and 
the  log  r  curve  for  points  in  the  corona  region  determines  rQ,  as  is  shown  in 
Figure  7.  Further,  it  is  seen  each  time  that  at  r  =  0.5  m,  which  is  the  loca¬ 
tion  of  the  grounded  electrode,  the  voltage  V(r)  is  not  zero,  but  has  some 
finite  value,  V$.  This  is  because  there  was  a  200  a  resistor  attached  between 
this  electrode  and  pulser  ground. 


2.  THE  CYLINDRICALLY  SYMMETRIC  CORONA 

Figure  9  shows  a  diagram  which  will  serve  to  illustrate  both  the 
cylindrically  and  spherically  symmetric  corona.  It  is  assumed  that  the  inner 
electrode  of  radius  a  is  connected  to  a  high  voltage  pulser  which  has  negli¬ 
gible  series  inductance  and  resistance.  The  corona  region,  for  which 
a  1  r  i.  ro>  is  characterized  by  a  resistivity  pi  and  permittivity  e2.  The 
ohmic  region  for  rQ  £  r  b  has  an  ohmic  resistivity  pQ  and  a  permittivity  e2. 
The  outer  conductor  of  radius  b  is  connected  to  pulser  ground  through  a 
resistor  R. 

It  is  a  simple  matter  to  apply  Gauss'  law  to  this  system  and  the  details 
will  be  omitted.  Note,  however,  the  following  points: 

a.  By  assumption,  the  electric  field  E  =  r  E  (r)  is  radial  and  pro¬ 
portional  to  r  1  in  both  regions.  This  leads  to  7  •  E  =  0.  Also,  the  resis¬ 
tivity  is  constant  in  each  region. 

b.  The  boundary  condition  at  rQ  is  that  the  normal  component  of 
the  current  density  is  continuous  across  the  boundary.  This  will  require  a 
charge  on  the  boundary  determined  by  the  permittivity  and  resistivity  of  both 
regions. 
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c.  If  both  the  electric  field  E(r)  and  the  voltage  V(r)  are 
expressed  as  a  function  of  the  voltage  on  the  inner  conductor  VQ,  the  permit¬ 
tivities  ex  and  e2  will  not  explicitly  appear  in  the  expressions  for  either  of 
these  quantities. 

d.  The  potential  at  r  =  b  is  Vs.  This  could  be  attributed  to  a 
resistor  between  this  electrode  and  ground  or  to  a  contact  potential  at  r  =  b. 

With  these  comments,  the  expressions  for  E(r)  and  V(r)  are: 
a  <  r  <  ro 


V(r)  =  [*n  L.  +  £l  in  +  V, 

0  L  ro  Po  r  J 

E(r)  £l  I 

D  p0  r 

r0  <  r  <  b 


(6) 

(7) 


Vn-V- 
V(r)  =_2_1 
D 


(8) 


E(r) 


Vvs  1 


where 


(9) 


D  = 


(10) 


We  may  obtain  a  number  of  useful  expressions  from  these  equations.  The 
critical  electric  field  (Eg)  may  be  obtained  from  Equation  9  by  setting  r  =  rQ: 
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E 


B  = 


Vvs 

D 


(ID 


If  a  is  the  angle  of  the  cylindrical  wedge  and  L  is  its  depth,,  then  the 
total  current  through  the  sample  I  is 


I 


ar0L-B 

'o 


aL  VVS 

*0  o 


Thus,  from  Equation  8  and  12 


(12) 


V ( r )  =  —  I  £n  —  +  V- 
aL  r  5 


(13) 


and 


av 

a(log  r) 


=  -  2.303 


(14) 


where  nn  10  =  2.303.  Thus,  if  3V/3  (log  r)  and  I  are  obtained  froin  the  exper¬ 
imental  data;  use  Equation  14  to  calculate  p0. 


From  Equations  6  and  8 


av 

3  log  r 


K  £l  a  <  r  <  rn 
Po  “  “ 

K  r0  <  r  <  b 


(15) 


where  K  =  -2.303  (VQ  -  Vs)/D.  Thus,  the  ratio  of  the  slopes  of  the  two  solid 
curves  in  Figure  7  will  give  pi/p0. 
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Finally,  it  is  of  interest  to  define 

x  =  r0/a  (16) 

and 

y  =  (V0  -  vs)/a  EB  (17) 

Then  rewrite  Equation  11 

y  =  -x  in  fix  +  £l  x  inx  (18) 

b  .  ^0 

We  have  plotted  y  versus  x  with  b/a  =  50  and  Pi/p0  as  the  parameter  in  Figure 
10.  For  small  pi/p0,  y  exhibits  a  maximum  at  some  x.  For  example,  for 
p\/p0  =  0,  the  maximum  will  occur  at  x  =  (b/a)e  1  =  18.4  which  corresponds  to 

rQ  =  b/e.  As  V0  increases  so  that  y  increases  with  Eg  assumed  to  be 

constant,  r0  will  increase.  If  VQ  is  so  large  than  rQ  or  x  reaches  the  point 
corresponding  to  the  maximum,  the  system  will  become  unstable  and  breakdown 
will  occur  all  the  way  to  r  =  b  (Ref.  8).  Thus,  the  corona  will  attach  to  the 
grounded  conductor. 

As  pi/p0  increases  from  zero,  the  maximum  moves  to  the  right  in  Figure  10. 
At  pj/Pq  =  0.203,  it  has  just  arrived  at  x  =  50,  or  r  =  b  with  y  =  39.7.  The 

system  is,  therefore,  unconditional ly  stable  for  Pi/p0  >  0.203  for  b/a  =  50. 

The  external  resistance  R  in  Figure  9  will  also  tend  to  stabilize  the 
corona.  It  is  shown  in  Reference  8  that  if  a  resistance  of  value 

Rc  =  p0/aL  (19) 

is  inserted  in  series  with  the  soil  sample,  the  corona  will  be  unconditionally 
stable.  For  a»/2  wedge  with  L  =  0.1  m  and  pQ  =  350  Ji-m,  Rc  =  2229  a. 

3.  THE  SPHERICALLY  SYMMETRIC  CORONA 

Figure  9  suffices  to  describe  this  system.  Assume  here  that  E  is  again 
radial  and  proportional  to  r  2.  Once  again,  7  •  E  =  0  and  px  is  a  constant. 
Points  2,  3,  and  4  in  the  previous  section  once  again  apply.  The  expressions 
for  V(r)  and  E(r)  are: 
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a  >  r  >  r0 


V(r).  =  V°~V,S- 
v  D' 


r/i . pi\jl .!+ £i  ii  + 

L\  p°/r°  b  p°  rJ 


E(r)  =  l2lll  BX  i_ 
O'  Po  r2 


rQ  >  r  >  b 


Hi-i) 


Vp-Vs  i_ 

D1  p2 


where 


>'  *  (\  -  - 
\  po)  ro 


I  +  £l  I 

b  p0  a 


From  Equation  22 


V  -V  1 
0  vs  1 


i —  1  +  £1  i  r. 


b  p0  a 


The  total  current  I  through  the  sample  is 


I  .  4»r02  fl 


so  that  the  total  resistance  from  the  high  voltage  electrode  to  r  =  b  is: 
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r  =  o' 

I  4it 

For  the  case,  b  ♦  «,  Pi/p0  0,  D‘  =  rQ  1 ,  so  that 

R-J®. 

4ffr0 

If  we  define 


(27) 


(28) 


x  =  r0/a 
and 

y  -  (V0  -  Vs)/aEB 
Then  Equation  24  becomes 


The  quantity  y  has  been  plotted  versus  x  with  b/a  »  50  and  pi/p0  as  the  param¬ 
eter  in  Figure  11.  The  system  has  an  instability  for  pi/p0  =  0  at  rQ  =  b/2. 
However,  it  quickly  stabilizes  at  small  values  of  Pi/p0*  For  example,  the 
system  is  stable  at  pi/p0  =  0.01. 

For  the  case  b  ■*•  •,  Equation  28  becomes 

y "  (‘  ‘  %)x  +  1,2  <30) 
and  the  system  is  unconditionally  stable. 

4.  EVOLUTION  OF  THE  ELECTROSTATIC  PHASE 

Soil  will  break  down  at  a  metal -soil  interface  when  the  electric  field  is 
about  1  MV/m.  For  example,  in  Figure  12  for  Exp.  A  where  the  electric  field 
was  initially  1.8  MV/m  (Fig.  B.ll,  Appendix  B) ,  it  can  be  seen  that  the  soil 
remains  ohmic  for  about  10  us,  and  then  breaks  down  characterized  by  the 
current  starting  to  rise.  Call  the  delay  until  breakdown  occurs,  t^ ,  as 
shown  in  Figure  12. 


ity  curves  for  a  spherical  soil  sample  that  is  symmetrically  broken 
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ure  12.  Current  versus  time  for  shot  012  (cylindrical  wedge  soil  sample).  The  soil 
remains  nearly  ohmic  until  time  t,.  after  the  application  of  Vrt. 


Breakdown  is  characterized  by  the  formation  of  the  corona,  a  subsequent 
lowering  of  the  sample  impedance,  and  an  increase  in  the  total  sample  current 
I  as  shown  in  Figure  13  with  the  accompanying  voltage  characteristic  on  the 
center  conductor  being  illustrated  in  Figure  14.  If  the  voltage  is  not  large 
enough  to  push  rQ  either  to  a  radius  where  instability  will  occur  or  to  the 
outer  conductor,  the  current  will  reach  a  maximum  at  a  time  tm  as  shown  in 
Figure  13  and  then  will  decline  to  zero.  If,  on  the  other  hand,  as  insta¬ 
bility  radius  is  reached  or  the  corona  reaches  r  =  b,  then  breakdown  to  the 
outer  conductor  at  r  =  b  occurs,  at  some  time  tB,  and  the  current  I  will 
increase  suddenly  at  tg  as  shown  in  Figure  15. 

In  Figure  12  where  p0  =  1500  O-m,  for  t  <  t^ ,  V0/I  =  16.1  kn  whereas  the 
calculated  ohmic  resistance  RQ  given  by 

D  P0*n  (b/a) 

Rq  =  -  (31) 

is  14.7  kn.  The  measured  resistance  is  characteristically  larger  than  the 
bulk  low  voltage  resistance  apparently  due  to  "a  small  increase  in  effective 
resistivity  with  increase  in  current  density,  a  small  increase  in  impedance 
near  the  small  diameter  center  electrode,  and  a  small  increase  with  time  of 
the  contact  impedance"  at  the  center  electrode  (Ref.  2).  This  subject  has 
been  treated  exhaustively.*  Also  the  relationship  between  the  applied 
electric  field  at  the  center  conductor  E0  and  the  time  of  initiation  of  break 
down,  t.j  has  been  shown  experimentally.  In  Figure  16,  the  dependence  of 
tj  with  E0  is  shown  where  the  parameter  is  geometry  and  in  Figure  17,  the 
parameter  is  soil  type.  In  all  cases,  t^  decreases  as  E0  increases.  Typi¬ 
cally,  t.,-  has  dropped  to  about  1  ys  if  E0  is  greater  than  2  MV/m. 

Only  the  stable  case  where  breakdown  does  not  reach  r  =  b  will  now  be 
considered,  and  the  current  reaches  a  maximum  as  in  Figure  13;  note  that  in 
Exp.  A,  ranges  from  100  to  100  ys.  There  is  a  rough  correlation  between 
tm  and  the  peak  applied  voltage;  tm  increases  with  increasing  V0  peak. 


♦Buried  Conductor  Studies,  AFWL-NTA-TN-83-0308.  Electrical  Breakdown  Charac¬ 
teristics  in  0.8  and  1.0  Meter  MX  Soil  Samples,  JAYC0R,  San  Diego,  Calif,  1981; 
and  Theoretical  Note-318,  1982. 
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jre  13.  Current  versus  time  in  shot  038  (cylindrical  wedge  soil  sample).  The  current 
rises  to  a  maximum  at  time  t_. 


038  09 
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Figure  15.  Curren 
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It  has  been  pointed  out  in  paragraph  2  that  given  a  plot  of  V(r)  versus 
log  r,  as  shown  in  Figure  7,  rQ  may  be  obtained  from  the  intersection  of  the 
log  r  approximations  to  the  ohmic  region  and  the  corona  region.  From  the 
slopes  of  these  two  curves,  Pi/p0  may  be  found  from  Equation  15  and  then 
Eg  may  be  calculated  from  Equation  11.  There  were  a  large  number  of  shots 
obtained  in  the  Exp.  A  in  a  350  Q-m  sand  mix  which  allows  these  quantities 
to  be  investigated  in  detail. 

Use  shot  399  as  an  example.  Plots  of  V/I  versus  log  r  for  this  shot  are 
shown  in  Figure  2,  and  V  versus  log  r  in  Figures  5  to  8.  Maximum  current  for 
this  shot  is  15.25  A  occurring  at  200  us.  From  these  figures,  as  time 
increases  rQ  increases  and  pi/p0  decreases.  At  30  us:  r0  =  0.127  m, 

Pi/p0  =  0.334;  at  60  us:  rQ  =  0.135  m,  pi/p0  =  0.28  at  100  us,  rQ  =  0.141, 
Pi/p0  =  0.23;  and  at  200  us,  rQ  =  0.16  m,  pi/p0  =  0.13.  For  time  greater 
than  200  us,  the  dynamic  region  considered  in  Section  V,  is  reached. 

In  the  other  shots  which  have  been  investigated,  rQ,  pi/p0,  and  Eg  at 
current  maximum  have  been  found.  For  all  these  shots  (12  in  all,  see  Figures 

■m. 

B.l  through  B.10  and  B. 12  through  B.31)  Eg  ranges  from  about  100  to  250  kV/m. 
It  is  interesting  that  this  range  coincides  with  that  found  in  Reference  6. 

The  correlation  between  Eg  and  p\/pQ  with  various  parameters  such  as  the  total 
current  I,  the  current  density  and  the  peak  VQ  have  been  investigated.  It  was 
found  that  these  quantities  correlate  best  with  V0  peak.  This  may  be  justi¬ 
fied  by  noting  that  the  breakdown  behavior  of  the  soil  may  be  determined  by 
the  initial  value  of  the  electric  field  at  the  center  conductor.  Note  that 
this  electric  field  E0  is  related  to  V0  peak  by 


k)/(a  “  9 


In  Figures  18  and  19,  Eg  and  p\/p0  versus  E0  have  been  plotted.  This  data 
has  been  processed  by  the  least  squares  linear  regression  method,  which  shows 
tnat  the  following  linear  relationships  best  fit  the  data: 

Ed  =  54.5  +  87.6  En,  R2  =  0.69  (33) 


o,/p„  =  0.97  -  0.48  E„,  R  =  0.87 


Figure  18.  Linear  regression  of  critical  electric  field  at  current  maximum  with  peak 
electric  field  at  center  conductor  p„  =  350  fi-m  sand  mix. 
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If  Eq  is  in  MV/m,  Eg  will  have  the  units  of  kV/m;  similarly,  EQ  should  be  in 
MV/m  in  Equation  34.  The  coefficient  of  determination,  R2 ,  is  shown  for  each 
least  squares  linear  fit.  The  range  of  E0  for  which  this  data  is  valid  is 
0.5  >.  E0  >  2  MV/m. 

In  Exp.  A,  there  were  some  data  available  from  which  plots  of  Eg  and 
Pi/pQ  versus  E0  were  obtained,  for  soils  with  other  ohmic  resistivities  than 
350  fl-m  (Figs.  B.32  through  B.45,  Appendix  B).  These  are  shown  in  Figures  20 
through  25  for  soils  with  resistivities  of  200,  170  and  130  fl-m.  These  data 
should  be  used  with  caution  because  generally  no  more  than  four  points  were 
available  to  determine  each  curve,  and  some  of  the  data  points  did  not 
correspond  to  tm  when  I  as  a  maximum,  but  were  obtained  at  an  earlier  time. 
With  this  qual ification ,  these  curves  tend  to  show  that  Eg  increases  as  the 
ohmic  resistivity  increases. 

5.  CALCULATION  OF  MAXIMUM  CURRENT 

One  of  the  real  achievements  of  the  electrostatic  model  is  its  ability  to 
predict  the  maximum  value  of  the  current  through  the  sample.  The  maximum 
current,  IM,  which  is  shown  in  Figure  13,  has  been  measured  for  a  number  of 
shots  for  the  350  fl-m  sand  mix. 

The  general  expression  for  current  through  the  sample  is  given  by 
<xL  rn 

i=-_^-eb  {35) 

If  r0  and  Eg  at  I  =  IM  are  used.  Equation  35  yields  IM.  But,  Eg  for  IM  is 
given  in  Figure  18  and  rQ  was  also  tabulated  for  this  current. 

In  Figure  26,  the  measured  value  of  IM  versus  V0  has  been  plotted  at  tm, 
i.e.,  the  center  conductor  voltage  at  the  time  when  the  current  is  maximum. 

Of  the  12  shots  for  pQ  =  350  fl-m,  4  were  for  a  sample  with  L  =  0.1  m,  while  8 
were  for  L  =  0.3  m.  The  experimentally  measured  values  of  IM  are  shown  by  the 
solid  dots  in  Figure  26. 

Alternately,  Equation  35  has  been  used *to  predict  IM  where  we  have  used 
the  actual  rQ  and  Eg  obtained  for  a  shot.  This  prediction  of  IM  is  shown  by 
the  open  circles  in  Figure  26.  The  agreement  is  usually  very  close.  The  per¬ 
centage  differences  are  tabulated  in  Table  1. 


Linear  regression  of  critical  electric  field  at  current  maximum 
electric  field  at  center  conductor.  =  200  fl-m  sand  mix. 


Linear  regression  of  (>J{>0  at  current  maximum  with  peak  electric  field  at 
center  conductor.  p0  =  200  n-m  sand  mix. 


regression  of  critical  electric  field  at  current  maximum  with  peak  electric 
t  center  conductor.  p„  =  170  fi-m  sand  mix  and  Soil-B. 


regression  of  critical  electric  field  at  current  maximum  with  peak  electric 
t  center  conductor.  =  130  fi-m  Soil-B. 


Linear  regression  of  pjf>0  at  current  maximum  with  peak  electric  field  at  center 
conductor.  =  130  n-m  Soil-B. 


t.ts.r?* 


current  versus  voltage  at  center  conductor.  Dots:  measured 
les:  calculated. 
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The  most  important  conclusion  to  be  drawn  from  the  generally  remarkable 
correspondence  between  the  measured  and  calculated  1^  is  that  the  corona  must 
be  very  azimuthally  uniform.  The  radial  probes  which  measure  V,  from  which 
rQ  and  Eg  are  determined,  lie  at  an  angle  of  45  deg  as  shown  in  Figure  1. 

Thus,  rQ  and  Eg  are  really  only  valid  for  this  angle.  If,  for  example,  they 
were  quite  different  at  some  other  angle,  then  the  current  density  would  be 
azimuthally  asymmetric  and  the  total  current  could  not  be  predicted  by 
Equation  35.  That  may  explain  the  discrepancies  (15.5  and  21.1  percent)  in 
shots  424  and  426.  The  fact  that  there  is  so  little  percent  difference  leads 
to  the  conclusion  that  the  corona  is  azimuthally  symmetric  and  rQ  and  Eg  do 
not  vary  with  angle.  The  fact  that  electrical  phenomena  within  the  90  deg 
cylindrical  wedge  sample  must  have  exact  images  in  the  wedges  constituting  the 
other  three  quadrants  of  the  circle,  as  well  as  in  the  wedges  stacked  on  top 
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6.  SPHERICAL  CASE 

Referring  to  Figure  9,  consider  a  spherical  geometry  in  which  the  spheri¬ 
cal  inner  electrode  has  a  radius  a  and  the  outer  spherical  shell  will  be 
removed  to  infinity,  so  that  b  +  Then,  Equation  25  may  be  written  as: 


(1  -  pl/Po)  i.  +  (£i)  i 


For  the  high  voltages  used  in  Exp.  C  (point-plane  experiments),  e.g.,  about 
200  kV,  the  estimate  of  p\/pQ  indicates  that  it  will  be  less  than  0.1.  There¬ 
fore,  if  Pi/p0  is  set  equal  to  zero  in  Equation  36  the  simple  result  is: 

ro  ■  VEB  <37> 

By  using  Eg  =  200  kV,  and  VQ  =  200  kV,  Equation  37  indicates  that  rQ  is  about 

1  m. 

As  an  indication  of  the  validity  of  these  calculations,  it  should  be 
noted  that  Exp.  C*  has  not  yielded  a  direct  measurement  of  the  radius  of  the 
spherical  corona  region.  Rather,  what  has  been  observed  is  that  while  the 
ohmic  resistance  from  the  high  voltage  electrode  to  infinite  ground  is  about 
200  a,  when  200  kV  is  applied,  the  resistance  initially  at  t  =  0  is  about  30  a 
and  then  declines  to  about  4  0.  Subsequently,  it  drops  to  1  0  where  it  is 
believed  that  a  streamer  is  propagating.  If  a  spherical  corona  region  which 
is  growing  in  radius  is  postulated,  one  may  thus  account  for  the  decreased 
resistance. 

The  total  current  flowing  in  the  soil  may  be  calculated  by  reference  to 
Figure  9.  At  r0,  it  is: 


Er 

I  =  _£  .  2 

o  4ir  rn 


♦Private  Communication,  Charles  Mallon,  JAYCOR,  San  Diego,  Calif. 
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Making  use  of  Equation  37,  the  resistance  R  to  infinite  ground  is: 


R  -fofB 

I  4it  v0 


(38) 


For  V0  having  its  peak  value,  200  kV,  and  using  the  measured  ohmic  soil  resis¬ 
tivity  at  the  experimental  site  for  Experiment  C, 


pQ  =  50  Q-m 

we  obtain: 

R  =  2.0  x  10’2  Eg(kV/m)  a 

A  resistance  of  4  si  could  be  accounted  for  if  Eg  =  200  kV/m  corresponding  to 
rQ  =  i  m. 
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V.  DYNAMIC  THEORY  OF  THE  CORONA  REGION 

In  Section  IV,  the  formation  of  the  corona  region  about  a  high  voltage 
electrode  was  discussed.  The  model  is  a  simple  one  in  which  there  are  two 
regions:  one  adjacent  to  the  high  voltage  electrode  called  the  corona  region 
and  the  second  called  the  ohmic  region.  The  corona  region  has  radius  r0  and 
resistivity  px ,  while  the  ohmic  region  exists  for  r  rQ,  and  has  the  ohmic 
resistivity  pQ.  The  corona  region  is  in  some  fashion  broken  down  electrically 

1.  THE  DYNAMIC  MODEL 

If  the  relaxation  phase  of  the  voltage  pulse  is  considered;  i.e.,  the 
time  after  the  current  has  reached  its  peak  value  and  is  declining  to  zero, 
then  a  third  region  can  be  identified  with  a  resistivity  p2,  which  is  depen¬ 
dent  on  a  current  density  J,  and  is  intermediate  between  px  and  pQ  (Ref.  6). 

By  introducing  this  third  region,  Liew,  et.  al.,  could  phenomenologically  pre¬ 
dict  their  experimental  results  quite  accurately.  They  referred  to  this  model 
as  a  dynamic  model  of  soil  breakdown. 

In  a  manner  similar  to  that  of  Reference  6,  the  following  model  of  corona 
relaxation  is  proposed.  As  the  voltage  rises  at  the  high  voltage  electrode, 
assume  that  the  soil  breaks  down  as  described  in  Section  IV;  i.e.,  the  soil 
breaks  down  out  to  a  maximum  radius  rQ  determined  by  Eg,  the  peak  applied 
voltage,  pL  for  r  _<  rQ,  and  pQ  for  r  >  rQ.  It  should  be  noted  that  at  rQ  in 
the  ohmic  region,  there  is  a  current  density 

jc  =  Eb/»o  <39> 

which  has  been  called  the  critical  current  density. 

As  the  applied  voltage  and  current  relaxes  to  zero,  a  new  region  forms. 

It  is  defines  by  r;  _<  r  £  rQ,  where  rQ  is  determined  from  the  electrostatic 
phase--call  this  region  2.  The  radius  r^  is  defined  such  that  at  rx ,  the 
current  density  is  Jc.  The  resistivity  in  region  2,  p2 ,  depends  on  J  and, 
therefore,  on  r.  Since  p2,  does  not  equal  pQ  the  resistivity  of  this  region 
undergoes  hysteresis— cal  1  this  the  relaxation  phase. 
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During  the  relaxation  phase,  we  can  thus  identify  three  soil  regions 
which  are  shown  in  Figure  27.  (This  figure  could  represent  either  a  cylindri¬ 
cal  or  spherical  system.): 

Region  1,  corona  region,  a  <  r  <  r^ 
p  =  plt  a  constant;  J  >  Jc. 

Region  2,  hysteresis  region,  rx  _<  r  rQ: 
p  =  p2(r);  J  <  Jc  with  J  »  Jc  at  r  »  r1( 

Region  3,  ohmic  region:  rQ  £  r  £  b. 
p  =  pQ,  a  constant. 

2.  RELAXATION  THEORY 

Consider  a  cylindrical  system,  the  calculations  shall  be  compared  with 
the  results  of  shot  399,  of  Exp.  A. 

For  region  2,  use 

p2  »  pi  +  (p0  -  pi)d  -  J(0/Jc)2  *  d°) 

which  is  approximately  the  expression  that  was  empirically  determined  to  give 
the  best  agreement  between  theory  and  experiment  (Ref.  6). 

Note  that  at  rx ,  p2  =  pi  and  the  resistivity  is  continuous  across  the 
boundary  between  regions  1  and  2.  From  Equation  1,  Ep  will  be  continuous 
(since  Jp  must  be  continuous),  and  there  will  be  no  charge  on  this  boundary. 

At  r  =  rQ,  however,  p2  *  p0,  so  that  the  resistivity  and  Er  is  discontinuous 
across  the  bondary  between  regions  2  and  3.  There  must  be  charge  on  this 
boundary. 

Let  V0(t)  be  applied  voltage  at  the  center  electrode  of  radius  a.  Then  a 
total  current  I(t)  will  flow,  and 

J(r,t)  =  I(t)ort  (41) 

where  a  is  the  wedge  angle  and  %  is  the  length  of  the  cylindrical  sample.  For 
Figure  27,  a  =  2*  rad. 
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Cylindrical  or  spherical  system  for  dynamic  mi 
corona;  region  2:  hysteresis;  region  3:  ohmic 
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The  resistance  of  the  system  R(r)  which  is  defined  as  the  resistance 
between  radius  r  and  the  grounded  electrode  at  r  =  b  shall  be  calculated.  By 
this  definition,  R(a)  would  be  a  total  resistance  such  that  I  =  VQ/R(a). 

Consider  two  closely  spaced  cylindrical  shells  with  the  inner  shell  at 
radius  r  and  the  separation  of  the  shells  being  dr.  Then  the  resistance  of 
the  region  between  the  shells  is: 


dR  =  p 


R(r)  = 


=  f  p  4L 


It  is  not  difficult  to  obtain 


Region  3,  r0  £  r  <_  b 


(s) 


R(r)  =  in 


Region  2,  rx  _<  r  _<  rc 


(S) (i.y 

‘if -*)'■(*-$ 
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Region  1,  a  <  r  £  rL 


(45) 


where 

rx  =  I/Jco  i.  (46) 

3.  RESULTS 

Shot  399  is  reproduced  in  Figure  2  where  the  ordinate  is  just  R(r)  in 
kilo-ohms  and  the  abscissa  is  log  r.  For  this  shot,  pQ  =  350  n-m,  a  =  tt/2  rad, 
z  =  0.1  tn,  a  =  0.01  m,  and  b  =  0.5  m.  Each  curve  shot  399  is  at  a  different 

time  with  curve  A  occurring  at  the  earliest  time  and  curve  F  at  the  latest 

time.  The  times  for  each  curve  are  shown  in  Figure  2;  also  shown  is  the 
current  I  that  was  measured  at  each  time. 

Under  the  assumption  that  Jc  is  a  constant,  maximum  radius  rQ  which  marks 
the  end  of  the  electrostatic  phase  would  occur  when  I  is  a  maximum.  This 
would  correspond  approximately  for  curve  B  at  200  us.  Measure  rQ  directly  by 
approximating  curve  B  by  two  straight  lines  which  are  seen  to  intersect  at 
r  =  0.16  m.  Take  this  value  of  r  to  be  rQ.  Using  the  techniques  described  in 
Section  IV,  para  4,  obtain  pi/p0  =  0.129  from  which  Eg  may  be  calculated  by 
the  method  described  to  be  231.5  kV/m.  Alternately,  since  I  was  measured  to 

be  15.25  A  for  curve  B,  at  rQ  =  0.16  m,  Jc  is  607  A/in’  ,  and  from  Equation  27, 

Eg  =  212.3  kV/m.  This  value  of  Eg  is  8  percent  lower  than  the  value  calcu¬ 
lated  using  these  techniques. 

Now  use  Jc  =  607  A/m2  and  rQ  =  0.16  m  to  predict  curves  C,  D,  E,  and  F. 
Adjust  pi/p0  as  may  be  necessary.  The  results  are  shown  in  Figures  28  through 
31  where  we  have  made  use  of  Equations  43  through  46.  The  value  of  Pi/p0 
required  are  as  follows: 
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CURVE  C 


Figure  28.  Dynamic  model  prediction  (dashed  line)  of  Curve  C  of  Figure  2,  p  /p 


prediction  (dashed  line)  of  Cu 


prediction  (dashed  line)  of  Curve  F  of  Figure  2 
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Curve 

I/A 

Pl/Po 

C 

15 

0.137 

D 

14 

0.121 

E 

10.75 

0.182 

F 

5 

0.409 

Thus,  a  good  agreement  may  be  obtained  between  the  dynamic  model  and  the 
experimental  data,  providing  that  pi/p0  is  allowed  to  increase  as  the  total 
current  decreases,  not  an  unexpected  result. 
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VI.  STREAMERS 

If  the  voltage  on  the  center  conductor  VQ  is  large  enough,  the  current 
through  the  sample  will  first  start  to  rise,  as  if  toward  a  maximum,  but  then 
will  suddenly  increase  as  is  shown  in  Figure  15.  This  sudden  increase  indi¬ 
cates  that  the  broken  down  portion  of  the  soil  sample  extends  all  the  way  to 
the  grounded  conductor,  or  that  it  has  attached  to  the  grounded  conductor. 

This  state  is  referred  to  as  attached  breakdown. 

Attached  breakdown  can  occur  for  the  following  reasons:  If  the  corona  is 
symmetric  at  its  boundary  as  is  the  case  in  Figure  26,  the  corona  will  attach 
if: 


a.  rQ  becomes  large  enough  so  that  it  equals  the  critical  stability 
radius;  if  the  system  is  stable,  or  the  critical  stability  radius  exceeds  b, 
the  corona  will  attach  if  VQ  is  so  large  that: 

b.  rQ  =  b. 

The  symmetric  corona  may  contain  streamers,  but  they  apparently  are  not  domi¬ 
nant  enough  to  significantly  distort  the  boundary,  so  that  the  boundary 
retains  the  symmetry  of  the  electrode  system.  Corona  attachment  is  discussed 
in  paragraph  1.  On  the  other  hand,  streamers  may  dominate  the  corona.  Then 
attachment  occurs  if: 

c.  The  streamer  propagates  as  described  in  Reference  7  which  charac¬ 
terizes  the  streamer  propagating  outward  with  a  velocity.  For  this  case, 
envision  the  streamer  pushing  the  corona  ahead  of  it  as  shown  in  Figure  4  so 
that  the  corona  attaches  first  followed  by  streamer  attachment.  Streamer 
attachment  is  discussed  in  paragraph  2. 

For  cases  a  and  b,  when  corona  attachment  occurs,  presume  that  the 
current  is  carried  more  or  less  uniformly  by  the  corona.  For  case  c,  when 
streamer  attachment  occurs,  the  current  is  presumed  to  be  carried  primarily 
by  the  streamer,  which  may  benefit  from  attachment  in  the  sense  that  its 
resistivity  may  decrease,  its  current  increase,  and  its  diameter  increase. 

This,  however,  is  speculation. 


* 
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As  seen  in  Figure  26  the  corona  in  Exp.  A  is  cylindrically  symmetric. 
Therefore,  corona  attachment  of  types  (a)  or  (b)  should  occur  in  these  experi¬ 
ments.  On  the  other  hand,  in  paragraph  2,  Exp.  B,  streamers  exist  after 
attachment  has  occurred.  Thus,  in  these  experiments,  the  attachment  should  be 
of  type  (c),  streamer  attachment. 

To  strengthen  these  conclusions,  refer  to  Figure  32  which  is  a  plot  of 
attachment  time  tB  versus  VQ  peak,  the  peak  voltage  on  the  center  conductor. 
These  data  were  obtained  from  both  Exp.  A  and  B.  The  x  and  open  circles  o  are 
data  from  Exp.  A.  It  is  seen  that,  at  least,  for  times  tg  less  than  1  ms 
attachment  does  not  occur  for  VQ  below  about  65  kV. 

These  data  were  obtained  from  a  system  in  which  a  =  0.01  m  and  b  =  0.5  m. 
If  the  data  for  Eg  and  Pi/p0  for  pQ  =  350  Q-m  in  Figures  18  and  19  are  used 
and  Equation  18,  it  is  possible  to  show  that  the  system  is  stable  up  to  V0  = 
65.7  kV.  At  that  voltage,  Eg  =  190  kV/m  and  pi/p0  =  0.17  leading  to  an  insta¬ 
bility  radius  of  0.4  m.  The  system  continues  to  be  unstable  for  VQ  larger 
than  65.7  kV.  This  result  is  in  accord  with  the  observation  that  in  Figure  32 
attachment  does  not  occur  for  Exp.  A  below  65  kV.  This  is  further  evidence 
that  corona  attachment  is  occurring  in  Exp.  A. 

On  the  other  hand,  data  from  Exp.  B  (represented  by  the  solid  dots)  shows 
a  correlation  between  VQ  and  tg  such  that  tg  clearly  decreases  as  VQ  increases, 
with  attachment  occurring  below  40  kV.  This  behavior  is  in  accord  with  cur¬ 
rent  theories  of  streamer  attachment.  For  example.  Smith*  predicts  a  phenom¬ 
enological  dependence  of  V0  and  tg  of  VQtgn  =  const  where  n  is  approximately 
3/8. 

1 .  CORONA  ATTACHMENT 

In  shot  302  (Exp.  A)  the  voltage  on  azimuthal ly  located  probes  is  mea¬ 
sured  (Figs.  B.46  through  B.54,  Appendix  B).  The  location  of  the  probes  are 

*  In  a  set  of  notes  which  have  not  been  published,  I.  Smith  predicts  a  phenom¬ 
enological  dependence  among  streamer  closure  distance  d(cm)  between  a  high 
voltage  electrode  and  a  ground  plane,  time  for  closure  t  (us),  and  voltage 
V  (kV)  of  the  high  voltage  electrode,  which  is 


v  t  3/8  d  1/6 

d 


12 


uu  j  uii  him  i 


+  360  x  1  mL,  0.003  m  Rod 
x  90  x  0.1  mL,  0.01  m  Rod 
o  90  x  0.3  mL,  0.01  m  Rod 


Figure  32.  Time  tg  for  attachment  to  occur,  x  and  o:  corona  attachment; 
streamer  attachment. 
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shown  in  Figure  1,  where  the  radius  of  the  probe  arc  is  0.144  m.  Also, 
b  =  0.5  m  and  a  =  0.01  m,  L  =  0.1  m.  A  flattened-out  view  of  this  arc  is 
shown  in  Figure  33  where  the  x  coordinate  measures  distance  along  the  arc,  and 
y  measures  distance  down  (Fig.  1)  into  the  soil  sample.  Since  L  =  0.1  m, 

0  <  y  <  0.1  m  while  x  ranges  from  zero  to  ir/2  x  0.144  =  0.226  m.  Figure  33 
shows  that  the  nine  probes  are  of  unequal  depth. 

In  Figure  34,  data  from  these  voltage  probes  are  presented  in  the  form  of 
a  resistance  Reff  which  is  defined  as: 

Reff  ■  (»0  -  »»)/!  <47> 

where  Vn  is  the  voltage  at  the  nth  probe  and  I  is  the  total  current  through 
the  sample.  Thus,  Reff  is  the  effective  resistance  between  the  center  con¬ 
ductor  and  probe  n  if  all  the  other  probes  were  at  the  same  potential  as  Vp. 
The  smaller  is  Reff»  the  closer  is  Vp  to  VQ.  The  voltage  Vn  versus  time  are 
given  in  Figures  B.12  through  B.20,  Appendix  B. 

The  current  I  through  the  sample  is  shown  in  Figure  15.  Attachment 
occurs  at  about  240  ps  and  that  from  300  us  to  about  700  ps,  the  current 
rapidly  decays.  At  time  greater  than  above  730  us,  the  soil  is  nearly  ohmic. 

From  Figure  34,  during  the  current  decay  (curves  for  300  ps  and  700  us) 
the  voltages  Vn  are  nearly  equal,  indicating  a  high  degree  of  azimuthal  sym¬ 
metry  of  the  voltage  at  0.144  m.  This  symmetry  strongly  supports  the  idea 
that  streamers  are  absent  and  that  corona  attachment  has  occurred. 

After  attachment,  the  resistance  of  the  sample,  R  =  V0/I,  is  shown  in 
Figure  35.  At  300  ps,  it  is  about  400  0,  and  then  increases  slowing  with 

time.  Assuming  corona  attachment,  the  resistivity  of  the  broken  down  soil  may 

be  calculated.  Assuming  that  it  is  uniform,  it  is  found  to  be  16.1  a-m.  The 
ohmic  resistivity  of  the  soil  sample  (Soil-B)  is  130  ft-m,  so  that  the  ratio  of 
these  resistivities  is  0.123.  From  Figure  25  this  ratio  is  typical  for  this 
type  of  soil. 

The  complete  picture  of  how  the  sample  resistance  R  varies  with  time  is 
shown  in  Figure  36  for  shot  403  (Figs.  8.55,  B.56,  Appendix  B).  The  corre¬ 
sponding  VQ  and  I  are  shown  in  Figures  B.21  and  B.22.  The  ohmic  resistance  of 
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ocation  of  azimuthal  probes  on  probe  arc  of  Figure  1  radius  of  probe  arc:  0.144 
measures  distance  along  arc,  y  into  soil  sample.  Streamer  location  also  shown. 


Figure  35.  Resistance  versus  time  for  shot  302  after  corona  attachment. 


the  sample  is  8717  ft.  As  the  current  starts  to  rise  after  ti ,  R  drops  rapidly 
to  below  2000  ft.  At  about  150  us,  attachment  occurs;  R  then  drops  below  400  a 
and  recovers  very  slowly.  Eventually  it  will  revert  to  its  ohmic  value. 

Returning  to  shot  302,  it  is  of  interest  to  examine  the  behavior  of  the 
sample  before  attachment.  In  Figure  15  the  current  just  prior  to  attachment 
(at  240  ps)  is  75  A.  This  current  may  be  predicted  from  Equation  13  if  we  use 
as  rQ  the  instability  radius.  The  peak  VQ  is  about  70  kV  or  E  =  1.8  MV/m. 

From  Figure  25  this  leads  to  pi/p0  =  0.1.  From  Figure  7  it  is  seen  that  the 
instability  radius  is  about  0.29  m.  Thus,  using  Eg  from  Figure  24,  we  cal¬ 
culated  from  Equation  13  that  I  =  91  A  compared  with  measured  75  A.  Consider¬ 
ing  the  error  in  Eg  due  to  the  few  points  available  and  other  factors,  the 
agreement  is  not  bad.  In  this  calculation,  a  symmetric  corona  was  assumed. 

In  Figure  34,  just  prior  to  attachment  at  200  us,  it  is  seen  that  Reff  is  azi- 
muthally  varying,  with  a  minimum  occurring  at  probe  2.  The  actual  voltages 
(at  200  us)  on  the  nine  probes  are  shown  in  Figure  37;  V0  at  200  us  is  65  kV.* 
This  suggests  that  there  may  be  a  streamer  near  probe  2  which  is  causing  its 
voltage  to  approach  VQ. 

A  model  for  the  corona  and  streamer  is  shown  in  Figure  38  just  prior  to 
attachment.  In  this  figure,  the  corona  boundary  has  arrived  at  the  insta¬ 
bility  radius,  0.29  m.  A  streamer  has  formed  at  the  center  conductor  and  has 
propagated  out  past  the  array  of  probes  and  is,  perhaps,  deforming  the  bound¬ 
ary.  Corona  attachment  subsequently  occurs.  It  is  interesting  that  within 
the  framework  of  this  model  that  the  streamer  clearly  disappears  as  is  shown 
in  Figure  34  for  300  us.  It  could  be  speculated  that  the  resistance  to  the 
outer  conductor  through  the  corona  is  lower  than  that  through  the  streamer. 
Thus,  the  corona  drains  current  from  the  streamer  and  the  streamer  disappears. 

It  is  possible  to  calculate  approximately  where  the  streamer  is  located 
and  estimate  its  radius.  An  extremely  simple  model  of  a  streamer  is  shown  in 
Figure  39.  The  streamer  is  modeled  as  a  cylinder  of  radius  aQ  with  a  resis¬ 
tivity  which  is  low  compared  to  the  surrounding  corona. 

♦This  number  was  calculated  from  the  data  in  Figures  34  and  37.  It  is  not 
clear  how  the  experimenters  measured  this  voltage. 


‘  V  V  /  - 


Figure  38.  Speculation  about  corona  boundary  and  streamer  location  in  shot  302 
just  prior  to  corona  attachment. 
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Let  Vs  be  the  potential  at  the  surface  of  the  streamer,  assumed  to  not 
differ  appreciably  from  the  center  conductor.  Then  if  V0 '  is  the  corona 
potential  in  the  absence  of  the  streamer,  the  potential  V(r)  in  the  vicinity 
of  the  streamer  may  be  approximated  as: 


V(r) 


W 


(48) 


where  r  is  the  cylindrical  coordinate  radial  distance  from  the  streamer  to  the 
point  where  V(r)  is  being  calculated,  and  bQ  is  some  reference  distance  where 
the  effect  of  the  streamer  is  no  longer  felt. 

Let  V(r)  be  evaluated  at  probes,  1,  2,  and  3,  which  appear  to  be  the 
closest  to  the  streamer.  Let 


Then  Equation  47  becomes 


(49) 


Vn  =  A  ln(b0/rn)  +  V0'  (50) 

where  rn  is  the  distance  from  the  streamer  to  the  nth  probe,  and  n  =  1,  2,  3. 

Call  (x,y)  the  location  of  the  center  of  the  streamer  where  it  cuts  the 
arc  shown  in  Figure  33.  Then  rn  =  rn(x,y).  Dividing  the  equation  for  and 

V3  by  the  equation  for  V2  (Eq.  50),  we  obtain  two  equations  which  do  not  con¬ 
tain  A: 


The  terms  Vi ,  V2 ,  and  V3  are  known  from  Figure  37,  and  we  estimate  V0' 
and  bQ.  Then  there  are  two  equations  and  two  unknowns,  x  and  y.  Recognizing 
that  rn  really  depends  on  xn'  shown  in  Figure  40  and  not  xn.  Equations  51  and 
52  have  been  solved  by  a  Newton-Raphson  technique.  We  have  used  VQ'  =  50.5  kV 
from  Figure  37  corresponding  to  V  at  a  remote  location  (probe  9)  and  for  bQ 
the  outer  radius  was  selected,  i.e.  bQ  =  0.5  -  0.144  =  0.356  m.  We  then  find 
that  x  =  0.26  m,  y  =  0.68  m.  The  location  of  the  streamer  is  shown  in  Figure 
33. 

These  values  of  x  and  y  will  guarantee  that  Vn,  n  =  1,  2,  3,  where  Vn  is 
given  by  Equation  50  will  have  the  values  shown  in  Figure  37.  Any  one  of 
these  equations  may  be  used  to  find  A:  A  =  3.2  kV.  Letting  Vs  s  VQ  =  65  kV, 

In  (bQ/a0)  =  4.53,  and  for  bQ  =  0.356  m,  aQ  =  0.0038  m  =  0.38  cm.  Thus,  the 
radius  of  the  streamer  is  predicted  to  be  0.38  cm. 

It  is  clear  that  there  are  many  sources  of  error  in  this  calculation, 
e.g.,  streamer  images  were  neglected,  and  the  value  of  bQ  is  unclear.  Thus, 
this  is  only  an  order  of  magnitude  calculation.  Note,  however,  that  the  pres¬ 
ence  of  streamers  can  be  easily  detected  by  probes  that,  as  in  Figure  33  are 
about  4  cm  apart. 

2.  STREAMER  ATTACHMENT 

As  an  example  of  streamer  attachment,  consider  shot  525  from  Exp.  B 
(Figs.  B.57  through  B.90,  Appendix  B).  In  this  shot,  the  parameters  of  the 
system  are  a  =  0.01  m,  b  =  1  m,  L  =  1  m,  and  pQ  =  300  fl-m  sand.  The  total 
current  I  through  the  sample  and  VQ  are  shown  versus  time  in  Figures  B.57  and 
B.58.  In  Figure  B.57  for  I,  an  error  occurred  in  recording  I  between  120  ys 
and  240  ys  so  that  the  peak  current  is  not  available.  The  plot  in  Figure  B.57 
is  truncated  during  this  time  interval.  Note  that  I  first  rises  slowly  until 
about  tB  =  80  ys  when  attachment  occurs.  The  current  rises  rapidly  to  its 
peak  value  (not  recorded)  and  then  decays  as  the  voltage  decays.  The  resis¬ 
tance  R  of  the  sample  is  shown  in  Figure  41.  It  drops  from  its  ohmic  value  of 
220  0  to  about  100  n  at  tB,  then  drops  rapidly  during  attachment  to  less  than 
10  O  and  then  rises  slowly  as  the  soil  recovers. 

Figure  42  is  a  plot  for  shot  525  similar  to  that  shown  in  Figure  34.  A 
set  of  azimuthal ly  located  voltage  probes  were  located  on  an  arc  with  a  6-in 
radius  at  the  top  insulating  plate  of  the  system.  The  probes  were  separated 
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Figure  41.  Resistance  versus  time  for  shot  525.  R  is  ohmic 
resistance.  1 
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by  11  deg;  there  were  19  probes  covering  198  deg,  approximately  a  semicircle. 
The  voltage  on  these  probes  versus  angle  at  50  us  have  been  plotted,  prior  to 
tg;  at  100  us  during  the  build-up  of  current  during  attachment;  and  at  200  us 
when  the  voltage  has  decayed  and  the  current  is  just  starting  to  decay.  This 
latter  time  would  correspond  to  the  time  300  us  in  Figure  34  when  it  was  seen 
that  the  voltage  was  uniform  with  angle.  For  shot  525,  however,  it  is  seen 
that  at  200  us,  the  voltage  is  not  at  all  uniform.  The  voltage  at  the  center 
conductor  V0  at  200  us  is  7.2  kv.  There  are  three  angles  at  which  the  probe 
voltage  is  close  to  7.2  kV:  -90  deg ,  -15  deg  and  65  deg.  This  is  interpreted 
to  mean  that  there  is  a  streamer  occurring  at  each  of  these  angles.  Since 
attachment  occurs  in  the  presence  of  streamers,  this  is  streamer  attachment. 

It  is  likely  that  the  streamer  carry  the  bulk  of  the  current. 
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VII.  SUMMARY 


1 .  SUMMARY 

The  important  points  made  in  the  report  are  listed  as  follows.  The 
electrode  configuration  is  cylindrical. 

a.  The  important  features  of  the  corona,  including  the  location  of  maxi¬ 
mum  rQ,  are  defined  by  about  30  us. 

b.  Sample  current  reaches  a  maximum  or  else  attachment  occurs  in  a  sym¬ 
metric  corona  in  about  100  -  200  ys. 

c.  The  corona  can  be  defined  by  two  parameters.  Eg  and  pi/p0.  These 
parameters  have  been  tabulated  here  at  maximum  current  and  find  that  both 
parameters  correlate  well  with  peak  VQ.  Note  that  maximum  rQ  which  occurs  at 
maximum  I  can  be  calculated  from  Eg  and  Pi/p0.  It  is  found  that  Eg  and  pi/p0 
depend  on  pQ.  We  have  tabulated  Eg  and  Pi/p0  versus  V0  peak  for  four  values 
of  p0. 

d.  For  the  cases  examined  here,  Pi/p0  “  0.8  -  0.1  and  Eg  was  100  -  200 
kV/m  at  maximum  I. 

e.  For  a  symmetric  corona,  maximum  I  could  be  accurately  predicted  from 
the  model,  requiring  only  that  peak  VQ  and  pQ  be  known. 

f.  The  relaxation  of  the  current  through  the  sample  as  the  voltage 
decayed  was  accurately  predicted  using  a  dynamic  model  in  which  a  third  region 
was  introduced.  It  was  required  that  Pi/p0  increase  slowly  from  its  value  at 
maximum  I.  Mo  general  data  on  this  increase  px/p0  are  presented. 

g.  Attachment  occurs  by  two  different  methods:  corona  and  streamer 
attachment. 

h.  Corona  attachment  is  associated  with  a  radius  of  instability  in  a 
cylindrical  system.  This  radius  can  be  calculated.  Also,  there  is  a  voltage 
which  can  be  calculated  below  which  the  system  is  corona  stable. 

i.  A  spherical  corona  is  inherently  stable. 
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j.  For  the  geometry  investigated  (b  =  0.5  m),  the  time  for  corona 
attachment  ranged  from  100  to  more  than  a  few  milliseconds  in  a  random  manner 
that  showed  no  correlation  with  pulser  voltage,  VQ. 

k.  Streamer  attachment  time  correlates  with  VQ. 

l.  The  location  of  streamers  can  be  predicted  by  voltage  probes  that  are 
spaced  4  cm  apart.  Probably  larger  probe  spacings  are  possible. 

m.  A  streamer  radius  is  of  the  order  of  0.4  cm. 

2.  DISCUSSION  OF  CORONA  ATTACHMENT 

This  is  a  phenomenon  which  has  not  been  recognized  previously.  Two  inde¬ 
pendent  pieces  of  evidence  which  support  its  existence  have  been  presented  in 
this  report.  The  conditions  under  which  this  type  of  attachment  exist  are 
more  restricted  than  for  streamer  attachment. 

a.  The  corona  must  have  a  high  degree  of  the  symmetry  of  the  initiating 
electrode  system.  This  was  the  case  in  the  ir/2  wedge.  It  has  been  theorized 
that  the  corona  and  streamer  images  in  the  other  three  quadrants  of  the  circle 
as  well  as  in  the  image  wedges  stacked  above  and  below  the  real  wedge,  contrib¬ 
uted  a  high  degree  of  symmetry  in  the  wedge.  Another  way  to  bring  about 
enhanced  symmetry  is  to  drive  the  center  electrode  at  a  voltage  so  high  that 
the  electric  field  at  the  surface  of  the  electrode  is  well  above  the  threshold 
value  for  metal -soil  interface  breakdown.  The  idea  here  is  that  if  the  corona 
really  consists  of  a  lot  of  small  streamers,  a  large  electric  field  will  cause 
many  of  these  to  form;  i.e.,  they  occur  with  a  high  density,  thus  keeping  the 
outer  surface  (corona  boundary)  uniform.  If,  on  the  other  hand,  only  a  few 
little  streamers  form,  the  boundary  is  poorly  defined,  and  there  is  no 
shielding  of  one  little  streamer  from  the  next.  Thus,  one  little  streamer 
will  start  to  propagate  as  described  in  Reference  7,  and  it  may  evolve  into  a 
large  streamer. 

b.  The  grounded  electrode  must  have  an  area  large  enough  so  that  it 
intercepts  a  significant  portion  of  the  corona  boundary  and  current  will  flow 
through  a  significant  portion  of  the  corona  volume. 
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Note,  that  if  a  symmetric  corona  arrives  at  a  grounded  electrode  which 
has  a  small  area  relative  to  the  corona  boundary,  it  can  be  speculated  that  a 
streamer  might  form  at  the  grounded  electrode  and  propagate  backwards  toward 
the  high  voltage  electrode.  Also,  note  that  in  most  of  the  experiments 
reported  to  date,  attachment  of  any  sort  has  not  occurred  over  distances  much 
larger  than  1  m.  Since  coronas  have  been  observed  to  have  radii  of  this  order 
of  magnitude,  they  should  not  be  ruled  out  as  the  main  mechanism  for  attach- 
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APPENDIX  A:  LINEAR  ELECTRICAL  CONDUCTIVITY  OF  SOIL* 


At  low  voltages,  soil  is  known  to  exhibit  electrical  conductivity 

5 

which  is  relatively  independent  of  frequency  below  10  Hz.  The  conductivity 
varies,  however,  from  soil  to  soil  and  varies  with  time  in  a  single  soil. 

This  suggests  that  the  conductivity  of  a  soil  is  a  complex  property  dependent 
on  a  number  of  variables.  The  purpose  of  this  Appendix  is  to  discover 
what  the  variables  are  and  to  attempt  to  define  them  quantitatively.  A 
literature  search  has  shown  that  the  variables  of  interest  include  the 
nature  of  the  soil  involved,  the  nature  and  amount  of  fluids  held  in  the 
soil,  and  the  temperature  of  the  system.  This  Appendix  will  discuss  these 
parameters  and  their  interrelationship  which  determines  the  conductivity 
of  a  given  soil. 

In  order  to  discuss  soil  conductivity,  we  must  first  define  what  is 
meant  by  'soil1.  A  simplistic  definition  of  soil  is  'weathered  rock’, 
or  rock  that  has  been  physically  broken  apart  and  chemically  and  biologically 
altered.  The  soil  is  a  complex  system  which  extends  from  the  surface  of  the 
earth  down  to  the  unaltered  rocks  of  the  earth's  curst,  and  consists  of 
minerals,  organic  material,  and  voids  between  grains  which  may  be  filled 
with  air  or  water.  Soil  has  a  diverse  spectrum  of  compositions  which  vary 
both  laterally  and  from  the  surface  to  extreme  depths.  In  general,  the  soil 
close  to  the  surface  has  a  high  proportion  of  organic  material,  such  as 
decaying  plant  remains,  organisms,  and  microorganisms,  while  the  soil  at  large 
depth  consists  almost  entirely  of  broken  up  bits  of  rock.  The  lower 
boundary  of  the  soil  is  a  point  of  debate  among  pedologists,  most  of  whom 

*  The  material  in  this  Appendix  was  prepared  by  Mr.  Rik  E.  Lantz,  a  graduate 
student  in  the  Department  of  Geosciences  at  the  University  of  Arizona.  It 
was  prepared  as  an  individual  studies  project  under  the  direction  of  one  of 
the  authors,  and  was  completed  in  December,  1982. 


claim  that  a  soil  must  have  at  least  some  organic  component  in  addition  to 
mineral  components  to  qualify  as  a  soil,  thus,  excluding  sand  and  gravel 
from  the  realm  of  soils.  Similarly,  pedologists  may  claim  that  lithified 
materials,  that  is,  materials  whose  component  grains  have  become  cemented 
together  by  precipitation  of  chemicals  carried  in  groundwater,  have  ceased 
to  be  soils  and  become  rocks.  From  an  electrical  standpoint,  however,  these 
are  philosophical  matters,  since  neither  organic  material  or  cements  have 
a  pronounced  effect  on  the  electrical  properties  of  the  materials  Involved. 

For  this  discussion,  the  soil  is  simply  a  veneer  on  the  earth's  surface 
between  the  actual  surface  and  unaltered  rock  at  large  depth. 

The  formation  of  soil  from  rock  is  a  complex  process,  involving  the 
physical  cracking  and  disintegration  of  the  rock  to  smaller  grains,  and 
chemical  activity  of  groundwater,  plants,  and  microorganisms  on  these 
grains.  Physical  breakdown  occurs  because  of  temperature  variations  at 
the  surface,  which  induce  stresses  between  grains  due  to  differing  coefficients 
of  thermal  expansion,  as  well  as  erosion  by  wind  and  water  and  a  myriad 
of  other  mechanical  processes.  Chemical  breakdown  occurs  because  most 
rocks  were  formed  at  pressures  and  temperatures  profoundly  different  from 
those  at  the  surface.  As  a  result,  the  mineral  phases  which  were  stable 
at  the  time  of  formation  are  not  stable  at  surface  temperatures  and  pressures, 
and  with  the  aid  of  water,  oxygen,  and  organic  activity,  they  break  down  to 
more  stable  phases,  notably  clays.  The  fact  that  weathering  begins  at  the 
surface  and  works  its  way  downward  has  marked  effects  on  soil  formation. 

As  mentioned  previously,  the  soil  near  the  surface  may  have  a  high  proportion 
of  organic  material.  Rainwater  which  falls  on  this  organic  matter  may 
become  acidic  and  begin  percolating  down  through  the  soil,  dissolving  certain 
solids  from  the  zone  immediately  below  the  organic  matter.  As  the  acidic 
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water  dissolves  more  and  more  solids,  its  pH  changes,  until  a  point  may  be 
reached  where  saturation  in  one  or  more  ions  occurs,  and  these  will  pre¬ 
cipitate. 

This  is  a  very  oversimplified  sketch  of  the  formation  of  soil  horizons, 
but  has  been  included  here  to  introduce  an  important  characteristic  of 
soil,  namely  that  it  is  exceptionally  heterogeneous.  Not  only  do  composi¬ 
tions  vary  from  one  location  to  another,  but  they  will  vary  with  depth  at 
a  single  location.  The  development  of  different  horizons  results  in 
horizontal  stratification  of  soil  and  may  juxtapose  two  types  of  soil  with 
contrasting  electrical  properties  in  a  relatively  short  vertical  distance. 

This  anisotropy  will  have  a  complex  effect  on  the  behavior  of  electric 
fields  in  the  soil  which  is  difficult  to  define  quantitatively.  The  follow¬ 
ing  discussion  has,  therefore,  been  limited  to  an  electrically  homogeneous, 
isotropic  soil,  and  the  relationships  discussed  will  apply  only  to  this 
case.  A  single  layer  within  a  stratified  soil  is  a  close  approximation 
to  such  a  case,  as  are  relatively  simple  porous  media,  such  as  loose  sand 
and  gravel. 

Soil  is  made  up  of  three  components:  mineral  material,  organic  material, 

and  voids  which  may  or  may  not  be  filled  with  fluid.  The  conductivities 

-1 2 

of  the  overwhelming  majority  of  minerals  is  very  low  (10  mhos/m  for 
quartz,  0.5  x  10  mhos/m  for  calcite,  10  mhos/m  for  muscovite),  as  is 
the  conductivity  of  dry  organic  matter.  As  a  result,  conductivity  in  the 
soil  is  dominated  by  the  third  component,  the  fluids  in  the  voids.  Pores 
in  the  soil  may  be  filled  to  varying  degrees  with  fluid,  usually  water 
with  some  dissolved  solids.  The  conductivity  of  wet  soil  is  much  greater 
than  the  conductivity  of  dry  soil.  Conductivity  of  a  slightly  damp  soil 

is  somewhere  in  between.  Put  simply,  the  more  water  available  to  conduct 
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electricity,  the  greater  the  conductivity  of  the  soil  as  a  whole.  Conductivity 
is  dependent  on  the  amount  of  space  available  to  be  filled  with  water  as 
well,  so  a  porous  garden  soil  will  have  a  greater  conductivity  than  a 
densely  compacted  chunk  of  clay.  The  conductivity  of  the  fluid  itself 
is  another  factor.  Finally,  conductivity  is  dependent  on  the  distance 
the  current  has  to  travel.  Soil  with  relatively  straight  water  filled 
pores  will  have  higher  conductivity  than  the  same  soil  with  contorted 
water  filled  pores.  The  conductivity  of  soil,  or  of  any  porous  medium, 
is  described  by  a  formula  derived  by  G.E.  Archie  in  1941,  which  takes  all 
of  the  above  factors  into  account.  This  relationship,  known  as  Archie's 
Law  [13],  is  as  follows: 
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where  a£  =  overall  conductivity  of  soil  with  fluid; 

a ^  *  conductivity  of  pore  filling  fluid; 

a  =  an  empirically  determined  constant  of  proportionality: 
0.5  <  a  <  2.5  ; 


<j>  =  porosity,  the  ratio  of  the  volume  of  the  voids  between  grains 
to  total  volume  of  soil  :0  <_  $  1 ; 

m  =  an  empirically  determined  constant  which  takes  into  account 
particle  shape  and  pore  geometry:  1.3  <_  m  <_  2 .5; 

s  *  fraction  of  pore  spaces  filled  with  fluid:0  <  s  <  1; 

n  *  an  empirically  determined  constant  n  ~  2. 

As  a  specific  case,  Archie's  Law  for  unconsolidated  sand  is  given  by: 
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where  a  *  0.65  [14] »  m  a  1 .5  [15],  and  n  =  2  [13]. 


Conduction  of  current  through  aqueous  solutions  is  achieved  by  ions 
which  carry  charge.  For  this  reason,  a  solution's  conductivity  is  governed 
by  the  concentration  and  electrical  affinity  of  these  ions.  Fig.  A.l  shows 
the  effect  of  concentration  of  various  ionic  solutions  on  resistivity. 

High  ion  concentrations  have  low  resistance  and  high  conductivity,  and  low 
concentrations  have  high  resistance  and  low  conductivity.  In  the  southwest, 
the  ion  content  of  groundwater  tends  to  be  very  high,  so  fluids  in  the 
soil  here  should  have  high  conductivity.  Solution  conductivity  is  also 
governed  by  temperature.  Low  temperatures  hinder  the  movement  of  ions  in 
solution  and,  thus,  the  solution's  capacity  to  carry  current.  A  rough 
rule  of  thumb  is  that  resistivity  changes  2%  per  degree  Centigrade  [16]. 

The  factor  <t  is  the  porosity  of  the  soil,  or  the  ratio  of  the  volume 
of  voids  (pore  spaces)  to  the  total  volume.  Porosity  ranges  from  0  to 
about  0.75,  and  depends  on  a  number  of  factors  including  grain  size,  range 
of  grain  sizes,  grain  shape,  and  packing.  Coarse  grained  soils  tend  to 
have  higher  porosities  than  fine  grained  soils.  A  wide  range  of  grain 
sizes  will  allow  very  small  grains  to  occupy  pore  spaces  and,  thus,  reduce 
porosity.  Platy  grains,  such  as  clays,  can  be  packed  together  more  closely 
than  nearly  spherical  grains,  such  as  rounded  quartz  sand  grains,  implying 
that  pure  clay  will  have  lower  porosity  than  sand.  Finally,  it  is  obvious 
that  a  compacted  soil  will  have  lower  porosity  than  an  uncompacted  soil. 

As  an  example,  Blatt  et  al .  [17]  quote  a  porosity  of  45%  for  freshly 
deposited  River  Sand.  Middleton,  et  al.  [18]  give  other  examples  of  porosity 
values  for  various  soil  types  which  are  shown  in  Table  Al . 
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Percentage  Solution 


Fig.  A.l  -  Diagram  showing  relationship  of  concentration  to  resistivity 
for  ionic  solutions  with  different  electrical  affinities.  (Fig.  1.2  from 
Tagg  [19]). 


TABLE  A1 


Kirvin  fine  sandy  loam  surface  29.7  76  53.4 

Taylor,  TX. 

Cecil  sandy  clay  loam  surface  44.0  69  52.0 

Statesville,  NC 

Shelby  silt  loam  surface  46.2  91  34.9 

Bethany,  MO 

Houston  black  clay  surface  61.0  89  40.2 

Temple,  TX 


The  factor  m  is  an  empirically  determined  constant  related  to  the  route 
the  current  travels  through  the  rock,  which  is  quantitatively  defined  as 
tortuosity.  Tortuosity  is  the  ratio  of  the  actual  length  of  a  path  through 
a  medium  to  the  linear  distance  between  the  starting  point  and  ending  point. 

The  more  tortuous  or  contorted  a  path,  the  greater  the  value  of  tortuosity. 

Due  to  the  complex  geometries  of  soil  pores  and  the  nature  of  electric 
fields  in  conductors,  current  can  flow  through  a  number  of  paths,  as  shown 
in  Fig.  A. 2.  The  tortuosity  in  this  case  is  simply  the  average  path  length 
the  current  flows  through,  divided  by  the  shortest  possible  path  (see  Fig. 

A.2).  Winsauer  et  al.  [14]  have  devised  a  simple  experimental  method 

* 

to  determine  tortuosity,  t,  directly,  and  have  found  that 

*  Winsauer  et  al .  assume  complete  saturation,  s  *  1  and,  in  this  case, 

o 

Archie's  Law  reduces  to  <j>  *  ~  .  They  state  in  equation  (9)  that 

1.67  ct  ae  a  ,  -7  . 

— - — *  ~  .  Substituting  <}>  for  ^  ,  the  equation  t  0/  =  $  "m  is  obtained. 
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Fig.  A. 2  -  Diagram  illustrating  the  flow  of  current  in  pores.  A  relatively 
small  amount  of  current  flows  through  the  straight  path  from  A  to  A'.  The 
average  of  all  of  the  current  paths  is  B  to  B ' ,  and  the  tortuosity  is 
BB'/AA'.  (Fig.  1  from  Winsauer  et  al.  [14]). 


Measuring  tortuosity  and  porosity  allows  empiric  determination  of  m. 

Tortuosity  is,  in  part,  a  result  of  particle  shape,  and  Jackson  et  al . 
[15]  have  attempted  to  evaluate  the  effect  of  particle  shape  on  m. 

Their  work  has  shown  that  particle  shape  and  m  are  directly  related,  and  that 
the  more  angular  the  particles,  the  greater  the  value  of  m.  The  result 
of  their  research  is  illustrated  in  Fig.  A. 3.  The  figure  shows  m  =  1.2 
for  commercially  produced  glass  spheres,  m  =  1 .85  for  extremely  angular 
shell  fragments  of  the  same  size,  and  intermediate  values  of  m  for  inter¬ 
mediate  grain  shapes.  Their  experimental  approach  to  empirical  determination 
of  m  was  to  completely  saturate  their  sands  with  fluids,  decrease  the 

porosity  by  successively  compressing  the  samples,  and  determine  m  as  the 

a 

slope  of  the  line  when  4>  was  plotted  against  ~  .  This  method  is  probably 

e 

more  accurate  than  the  method  derived  from  Jackson's  work  because  it  is  a 
more  direct  measurement. 

The  parameters  take  into  account  the  amount  of  fluid  actually  held 
in  the  pore  spaces.  The  parameter  is  a  saturation  factor  ranging  between 
0  and  1 ,  which  is  simply  a  ratio  of  the  amount  of  pore  space  filled  with 
fluid  to  the  total  pore  space.  The  parameter  is  yet  another  empirically 
determined  constant,  which  can  be  evaluated  by  determining  all  other  con¬ 
stants  and  varying  s.  There  are  some  physical  constraints  which  in  very 
dry  conditions  limit  the  values  of  s.  One  of  these  is  the  phenomenon 
of  adsorption,  which  results  from  an  electrical  imbalance  of  minerals 
in  the  soil.  Due  to  the  nature  of  crystal  growth,  the  surfaces  of  min¬ 
erals  often  have  unneutral ized  charges,  which  can  attract  one  or  more 
layers  of  dipolar  water  molecules  (see  Fig.  A. 4).  Even  when  the  mineral 


Fig.  A. 3  -  Diagram  showing  the  effect  of  particle  shape  on  m.  The 
particles  are  increasingly  angular  from  A  to  E.  The  linear  groups  of 
points  reflect  successive  decreases  in  porosity  for  each  sample.  Their 
notation  has  o  /a  -  FF,  and  <f>  =  n.  (Fig.  6  from  Jackson  et  al .  [15].) 


Fig.  A. 4  -  Diagram  showing  charge  distribution  at  the  surface  of  a 
mineral  crystal  and  two  adsorped  layers  of  water  molecules.  The  symbol 
=  +  represents  a  dipolar  water  molecule. 


Fig.  A. 5  -  Schematic  illustration  of  meniscus  contacts.  Water  is  in 
black,  grains  are  shaded,  and  air  filled  pore  space  is  white.  When  the 
radius  of  curvature  of  the  meniscus  gets  very  small,  surface  tersion  force 
are  high,  and  the  water  will  not  evaporate. 
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is  dried,  the  thin  lamina  of  water  remains  and  provides  a  path  for  current 
flow.  This  effect  is  especially  pronounced  in  clays.  Another  phenomenon 
which  keeps  soil  from  drying  completely  is  surface  tension.  At  low 
fluid  content,  water  Is  held  very  tightly  In  between  grains  in  what  are 
known  as  meniscus  contacts  due  to  surface  tension  (see  Fig. A.5).  A  certain 
minimum  value  exists  for  each  soil  where  surface  tension  forces  are  greater 
than  the  water's  natural  tendancy  to  evaporate,  and  the  water  held  in  these 
meniscus  contacts  will  remain  unless  the  soil  is  subjected  to  abnormal 
conditions.  The  result  of  these  two  physical  constraints  on  saturation 
is  that  each  soil  contains  a  minimum  amount  of  water  and,  therefore,  should 
have  a  minimum  conductivity.  This  minimum  value  can  be  significant,  as 
shown  in  Fig.  A. 6,  where  the  minimum  value  of  s  for  a  'clay  soil'  is  around 
10  per  cent. 

In  summary,  then,  it  has  been  shown  that  soil  conductivity  for 
frequencies  below  about  100,000  Hz  is  described  by  the  relatively  simple 
relationship,  Archie's  Law.  The  form  of  Archie's  Law  indicates  that  the 
solid  components  of  soil  are  more  or  less  electrically  Inert,  and  that 
the  parameters  which  govern  soil  conductivity  are  porosity,  amount  of 
pore  fluids,  ionic  content  of  pore  fluids,  and  grain  shape.  The  soil 
is  such  a  complex  system  that  the  law  must  be  modified  by  several  empirically 
determined  constants  which  allow  for  compositional  and  geometric  variations 
between  soils  in  order  to  describe  each  individual  case.  After  experimentally 


determining  these  constants,  an  equation  can  be  written  which  accurately 
expresses  the  conductivity  of  a  given  soil,  and  the  meaning  of  such  an 
equation  has  been  the  object  of  this  investigation. 
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Fig.  A. 6  -  Hysteresis  curves  showing  per  cent  moisture  vs.  pF,  a  measure 
of  the  free  energy  of  molecules  of  the  water  held  in  the  soil,  which  is 
proportional  to  the  amount  of  energy  needed  to  get  them  out.  Point  Y 
represents  'air  dryed'  soil,  and  shows  the  associated  minimum  water  content, 
a)  Loam  Soil;  b)  Sandy  Soil i  c)  Heavy  Clay  Soil.  (Fig.  4  from  Ingram,  [20]). 
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APPENDIX  B 
DATA  BASE 

Presented  in  this  Appendix  is  the  raw  data  from  the  experiments  upon 
which  many  of  the  calculations  in  this  report  are  based.  The  data  are  digita¬ 
lizations  of  Polaroid  photos  and  are  presented  as  plots  of  voltage  or  current 
versus  time. 

In  the  figure  captions,  the  first  number  is  the  shot  number  followed  by  a 
letter  designating  the  group  of  experiments  from  which  it  was  obtained.  The 
second  number  is  a,  the  radius  of  the  center  conductor;  the  third  is  b,  the 
radius  of  the  outer  conductor.  The  fourth  is  L,  the  sample  length,  the  fifth 
is  ohmic  resistivity  in  ohm-m.  All  distances  are  in  meters.  For  example, 

399 A,  0.01,  0.5,  0.1,  350  sand  mix  means  shot  399  of  Exp.  A,  with  a  =  0.01  m, 
b  =  6.5  m,  L  *  0.1  m,  and  pQ  =  350  n-m  sand  mix. 

Also,  V0  indicates  that  the  voltage  on  the  center  conductor  in  kilovolt 
is  being  plotted  versus  time  in  microsecond;  I  means  that  total  sample  current 
in  amperes  is  being  plotted.  Finally,  V  indicates  that  a  probe  voltage  is 
being  plotted;  and  r  is  the  radial  position  of  a  probe  in  meters;  e  is  azimuthal 
angle  of  probe  in  degrees;  and  bottom  or  top  indicates  that  a  probe  was  near 
the  bottom  on  top  of  the  soil  sample. 


Figure  B.3.  399A,  .01,  .1,  .5,  .1,  350  sand-mix;  V  at  radial  position  2;  r  =  .059  m. 


position 


Figure  B.6.  399A,  .01,  .1,  .5,  .1.  350  sand-mix;  V  at  radial  position  5;  r  =  .209  m. 


radial  position  6;  r  =  .259  m 
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Figure  B.8.  399A,  .01,  .1 ,  .5,  .1,  350  sand-mix;  V  at  radial  position  7;  r  =  .309  m. 


Figure  B.9.  399A,  .01  ,  .1  ,  .5,  .1  ,  350  sand-mix;  V  at  radial  position  8;  r  =  .359  m. 


399A,  .01  ,  .1  ,  .5,  .1  ,  350  sand-mix;  V  at  radial  position  9;  r  =  .409  m. 


1500  sand; 


Figure  B.12.  394A,  .01,  .5,  .1,  350  sand  mix; 


Figure  B.13.  394A,  .01,  .5,  .1,  350  sand  mix 


Figure  B.14.  401A,  .01  ,  .5,  .1  ,  350  sand  mix; 


120.  140.  100.  100.  200 


9999999999 

n  n  u  n  k 


Figure  B. 18.  420A,  .01,  .5,  .3,  350  sand  mix; 


350  sand  mix; 


Figure  B.26.  424A,  .01,  .5,  .3,  350  sand  mix; 


Figure  B.28.  426A,  .01,  .5,  .3,  350  sand  mix; 


Figure  B.29.  426A,  .01,  .5,  .3,  350  sand  mix; 


Figure  B.31.  427A,  .01,  .5,  .3,  350  sand  mix; 


Figure  B.32.  170A,  .01,  .5,  .1,  200  sand; 


Figure  B.33.  170A,  .01,  .5,  .1,  200  sand; 


Figure  B.35.  171  A,  .01  ,  .5,  .1  ,  200  sand; 


sand; 


200  sand; 


Figure  B.38.  268A,  .01,  .5,  .1,  170  Soil-B; 


Figure  B.40.  272A,  .01,  .5,  .1,  170  Soil-B; 


2  10  69KU 


Figure  B.41  .  272A,  .01, 


Figure  B.42.  281A,  .01,  .5,  .1,  130  Soil-B; 


Figure  B.44.  282A,  .01,  .5,  .1,  130  Soil-B; 


Figure  B.45.  282A,  .01,  .5,  .1,  130  Soil-B; 


Figure  B.46.  302A,  .01,  .5,  .1,  130  Soil-B;  V  at  azimuthal  position 


Figure  B.47.  302A,  .01,  .5,  .1,  130  Soil-B;  V  at  azimuthal  position 


Figure  B.48.  302A,  .01,  .5,  .1,  130  Soil-B;  V  at  azimuthal  position 


azimuthal  position 


Figure  B.bO.  302A,  .01,  .5,  .1,  130  Soil-B;  V  at  azimuthal  positi 


Figure  B.51.  302A,  .01,  .5,  .1,  130  Soil-B;  V  at  azimuthal  position 


Figure  B.52.  302A,  .01,  .5,  .1,  130  Soil-B;  V  at  azimuthal  position 


imuthal 


Figure  B.54.  302A,  .01,  .5,  .1,  130  Soil-B;  V  at  azimuthal  position 


Figure  B.55.  403A,  .01,  .5,  .1,  350  sand-mix; 


Figure  B.56.  403A,  .01  ,  .5,  .1  ,  350  sand-mix; 


B.59.  525B 


323 


Figure  B.60 


(bottom) 


at  r,0  =  .18  m,  0°  (bottom) 


.30  m,  0°  (bottom). 


(bottom) 


.15  m,  261  (top). 


525B,  .01  ,  1.0,  1.0,  300  sand-mix;  V  at  r,e  =  .15  m,  272  (top). 


Figure  B.68.  525B,  .01,  1.0,  1.0,  300  sand-mix;  V  at  r,e  =  .15  m,  283  (top) 


.15  m,  294  (top). 


525B,  .01  ,  1.0,  1.0,  300  sand-mix;  V  at  r,e  =  .15  m,  305  (top) 


525B,  .01  ,  1.0,  1.0,  300  sand-mix;  V  at  r,e  =  .15  m,  316  (top). 


Figure  B.72.  525B,  .01,  1.0,  1.0,  300  sand-mix;  V  at  r,e  =  .15  m,  327  (top) 


Figure  B.74.  525B,  .01  ,  1.0,  1.0,  300  sand-mix;  V  at  r,e  =  .15  m,  338  (top) 


.15m,  338  (top). 


Figure  B.76.  525B,  .01  ,  1.0,  1.0,  300  sand-mix;  V  at  r,e  =  .15  m,  349  (top) 


Figure  B.77. 


525  8A 


Figure  B.79. 


sand -mix;  V  at 


529  9B 


Figure  B.83.  525B,  .01  ,  1.0,  1.0,  300  sand-mix;  V  at  r,e  =  .15  m,  33  (top) 


Figure  B.84.  525B,  .01,  1.0,  1.0,  300  sand-mix;  V  at  r,e  =  .15  m,  44  (top). 


V  at  r,e  =  .15  m,  55  (top). 


525B. ,  .01  ,  1.0,  1.0,  300  sand-mix;  V  at  r,e  =  .15  m,  77  (top). 


Figure  B.88.  5258.,  .01  ,  1.0,  1.0,  300  sand-mix;  V  at  r,e  =  .15  m,  88  (top) 


300  sand-mix;  V  at  r,e  =  .15  m#99  (top). 


